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SUMMARY 


An  experimental  and  theoretical  research  progreun 
was  conducted  to  evaluate  the  performance  benefits  obtained 
from  the  vortex  flows  generated  in  the  pressure  field  of  a 
low  aspect  ratio  wing.  Wind  tunnel  studies  were  conducted 
with  an  unmodified  wing  panel  and  with  strakes,  snags,  sep¬ 
arator  plates  and  combinations  thereof  attached  to  the  panel. 
Additional  tests  were  conducted  to  determine  the  effect  of 
the  main  wing  wake  on  the  effectiveness  of  a  horizontal  tail 
surface.  The  results  of  the  test  program  showed  that  the 
performance  characteristics  of  the  unmodified  wing  panel  can 
be  improved  substantially  in  the  post-stall  region  with  the 
addition  of  leading  edge  devices  such  as  strakes.  The  re¬ 
sults  of  the  tests  with  tail  surface  indicated  that  while 
the  vortex  flows  developed  large  pitching  moments  on  the 
wing,  the  horizontal  tail  surface  placed  properly  in  the 
wake  was  able  to  counter  the  pitching  moment  and  trim  the 
configuration. 

Concurrently,  a  theoretical  analysis  was  developed 
based  on  a  modified  potential  lifting  surface  theory  and  a 
self-similar  solution  of  the  Navier-Stokes  equations  to  pre¬ 
dict  the  lift  and  pressures  on  the  wing  surface.  The  theory 
included  proper  accounting  of  the  potential  and  separated 
flows  on  the  wing  and  the  suction  forces  of  the  leading  edge 
vortices.  A  good  correlation  was  obtained  between  the  exper¬ 
imental  and  the  predicted  results  for  Navy's  F-4  type  wing 
and  a  reasonable  correlation  was  obtained  for  sharp  edged 
delta  wings. 
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FOREWORD 


The  work  described  in  the  technical  report  was  per¬ 
formed  by  the  RASA  Division  of  Systems  Research  Laboratories, 
Inc.,  for  the  Department  of  the  Navy,  Office  of  Naval  Research, 
Arlington,  Virginia,  under  Contract  Number  N00014-74-C-0091. 

The  research  program  was  underteiken  under  the  technical  cog¬ 
nizance  of  Dr.  Robert  E.  Whitehead  of  the  Vehicle  Technology 
Program  of  ONR. 
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1.  INTRODUCTION 


A  four  year  research  program  was  conducted  by  RASA 
Division  of  Systems  Research  Laboratories,  Inc.,  under  the 
sponsorship  of  the  Office  of  Naval  Research  to  study  the 
vortex  flow  interactions  with  the  lifting  surfaces  and  the 
possibility  of  applying  the  vortex  control  techniques  to 
improve  the  performance  and  maneuvering  capabilities  of 
present  generation  aircraft.  The  study  included  both  the 
experimental  measurements  of  performance  and  pressure  dis¬ 
tributions  on  the  lifting  surface  and  a  theoretical  analy¬ 
sis  of  the  flow  field. 

The  objective  of  the  research  progreim  in  the  first 
year  was  to  explore  the  feasibility  of  generating  a  signifi¬ 
cant  incremental  lift  through  wing/vortex  interactions.  A 
review  of  experimental  data  previously  obtained  on  the  F-4 
Phantom  wing  planform  showed  that  a  relatively  weak  snag 
vortex  generated  within  the  pressure  field  of  the  wing  can 
produce  a  large  beneficial  suction  peak  and  incremental 
lift.  Because  of  this  noted  apparent  benefit,  the  following 
investigations  were  undertaken  during  the  first  year  of 
research. 

1.  A  review  of  current  theoretical  capabilities 
in  the  field  of  vortex  control  and  wing/vortex  interactions. 

2.  Development  of  a  simplified  theoretical  model 
of  the  wing/vortex  interaction. 

3.  Study  of  "vortex  generator"  characteristics. 

4.  Comparison  of  representative  vortex  control 
configurations . 
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A  detailed  discussion  of  the  results  of  this  study  were 
presented  in  Reference  1. 

At  the  end  of  the  first  year’s  effort,  it  was  con¬ 
cluded  that  the  maximum  amount  of  incremental  lift  obtain¬ 
able  from  attached  vortex  flows  could  be  achieved  when  the 
vortex  was  close  to  the  lifting  surface  and  when  the  vortices 
were  formed  within  the  pressure  field  of  the  wing.  There¬ 
fore,  in  the  second  year,  an  experimental  program  was  under¬ 
taken,  in  part  to  demonstrate  and  to  verify  that  the  bene¬ 
ficial  incremental  lift  and  wing  vortex  interaction  mechanisms 
indicated  by  the  simplified  theory  utilized  in  Reference  1 
could  be  achieved. 

The  studies  were  conducted  using  an  F-4  type  of  wing 
planform  with  and  without  leading  edge  control  devices  in¬ 
cluding  a  strake,  a  snag,  and  combinations  thereof.  The 
pressures  and  the  total  forces  and  moments  on  the  wing  were 
measured.  Flow  visualization  studies  were  also  carried  out 
using  the  helium  bubble  technique.  The  data  obtained  was 
reduced  in  order  to: 

1.  amalyze  the  results  of  vortex  control  devices 
in  terms  of  improving  the  total  lift  and  drag  performance 
of  the  system  tested, 

2.  interrelate  the  improvements  in  lift  and  drag 
performance  with  wing/ vortex  interaction  phenomena  on  the 
basis  of  the  observed  spanwise  and  chordwise  pressure 
distributions , 

3.  obtain  a  more  basic  understanding  of  the 
parameters  which  control  the  positive  benefits  that  can  be 
achieved  through  the  interaction  of  a  vortex  with  a  lifting 
surface , 
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4.  determine  the  capability  of  the  available 
theoretical  model  in  relation  to  the  observed  experimental 
phenomena,  and 

5.  outline  the  theoretical  and/or  experimental 
research  requirements  to  further  develop  the  understanding 
of  the  use  of  vortex  flows  to  improve  the  performance  char¬ 
acteristics  of  the  aircraft. 

The  results  of  the  second  year's  efforts  and  a  detailed  dis¬ 
cussion  of  the  experimental  program  can  be  found  in  Ref¬ 
erence  2. 

These  results  in  general  demonstrated  that  the  most 
significant  effect  was  realized  by  the  addition  of  a  strake 
at  the  wing  root.  While  the  strake  vortex  by  itself  did  not 
generate  a  large  suction  lift,  it  strengthened  and  stabilized 
the  leading  edge  vortex  and  effected  a  significant  increase 
in  its  suction  lift.  The  addition  of  an  outboard  snag  re¬ 
sulted  in  a  further  beneficial  change  in  performance.  It 
was  also  concluded  that  the  effectiveness  of  the  root  strake 
had  reached  its  limit  and  additional  vortex  control  must  be 
sought  after  by  applying  devices  on  the  outboard  section  of 
the  wing.  Therefore,  an  additional  experimental  program  was 
undertedcen  in  the  third  year  to  investigate  whether  the  flow 
on  the  outboard  section  of  the  wing  could  be  stabilized.  In 
addition  to  that,  it  was  believed  that  the  beneficial  effects 
of  the  vortex  control  techniques  could  be  realized  at  lower 
angles  of  attack  by  generating  the  vortices  at  lower  angles 
of  attack  by  artificial  means.  It  was  also  concluded  from 
the  experimental  observations  of  the  second  year's  efforts 
that  the  simplified  theory  of  Reference  1  was  totally  in¬ 
adequate  to  handle  the  complex  flow  phenomenon.  Therefore, 
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additional  efforts  were  made  to  theoretically  represent  the 
complex  flow  field  over  the  lifting  surface  in  a  more  com¬ 
prehensive  and  concise  manner.  The  results  of  the  third 
year's  efforts  were  reported  in  Reference  3. 

The  test  results  obtained  at  this  stage  of  the 
research  effort  showed  that  while  positive  gains  in  lift 
and  drag  characteristics  could  be  achieved,  the  pitching 
moments  tended  to  become  unstable  with  increasing  angles 
at  higher  angles  of  attack.  Therefore,  an  experimental 
investigation  of  the  tail  effectiveness,  under  the  condi¬ 
tions  that  actually  existed  behind  a  lifting  surface  devel¬ 
oping  strong  attached  vortex  flows  was  undertaken  during 
the  fourth  year. 

The  results  of  a  doublet  lattice  method  developed 
during  the  third  year  to  predict  the  nonlinear  suction 
effects  showed  reasonable  comparisons  between  the  predicted 
and  observed  performance  and  pressure  loadings.  However, 
due  to  certain  inherent  assumptions  (such  as  neglecting 
the  thickness  and  the  curvature  effects)  made  in  the  analy¬ 
sis,  the  prediction  capabilities  of  that  analysis  were 
thought  to  be  rather  limited.  Therefore,  the  theoretical 
analysis  was  extended  to  include  the  effect  of  the  thick¬ 
ness  and  the  surface  slope.  Also,  a  better  model  of  the 
nonlinear  suction  effects  of  the  interacting  vortices  was 
incorporated  into  the  analysis.  The  results  of  the  efforts 
during  the  fourth  year  were  reported  in  Reference  4. 

In  the  fifth  year,  the  computer  program  written 
for  the  theoretical  analysis  was  streamlined  and  documented 


so  that  it  could  be  used  as  a  research  tool  to  analyze 
certain  other  planforms.  This  final  report  presents  a 
general  summary  of  the  entire  research  effort  and  the 
documentation  of  the  computer  program. 
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2.  TECHNICAL  DISCUSSION 


2.1  PRELIMINARY  INVESTIGATION  OF  THE  BASIC  MECHANISMS  OF 
WING/VORTEX  INTERACTION 

It  has  been  known,  for  several  years,  that  the 
interaction  between  the  wing  and  a  vortex  can  produce 
additional  lift  as  in  the  case  of  the  vortices  generated  by 
strakes,  close  coupled  canards,  and  snags.  The  vortices 
interacting  with  the  wing  can  be  free  vortices  created  by 
an  independent  upstream  aerodynamic  device  or  may  be  formed 
within  the  pressure  field  of  a  wing  surface  due  to  loading 
gradients.  The  chordwise  pressure  distributions  on  the 
upper  surface  of  a  symmetric  airfoil  in  the  vicinity  of  a 
free  vortex  are  shown  in  Figure  1  (Ref.  5)  to  illustrate 
the  incremental  suction  effect  from  a  free  vortex  generated 
by  an  upstream  lifting  surface.  As  shown  in  Figure  2  (Ref.  6)^ 
a  similar  incremental  suction  effect  is  also  observed  in  the 
case  of  a  snag  vortex  interacting  with  the  wing  surface  that 
generated  the  vortex. 

The  wing/vortex  interaction  basically  involves  two 
mechanisms.  The  induced  velocity  distribution  over  the  lift¬ 
ing  surface  produced  by  the  vortex  changes  the  effective 
angle  of  attack  and  the  loading  in  the  vicinity  of  the  vortex. 
The  change  in  the  loading  thus  obtained  is  such  that  the 
loading  increases  on  one  side  of  the  vorte.\  axis  (upwash 
side)  2md  decreases  on  the  other  side  (downwash  side)  of  the 
vortex  axis.  The  other  effect  is  the  nonlinear  suction  effect 
caused  by  the  low  pressure  within  the  rotational  field  of  the 
vortex.  Figure  3  illustrates  the  qualitative  nature  of  the 
vortex  suction  effect  for  a  pair  of  counter-rotating  vortices 
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on  the  suction  side  of  a  lifting  surface.  Based  on  simple 
physical  principles,  a  theoretical  model  was  developed  in 
Reference  1  to  predict  the  incremental  suction  peaks  pro¬ 
duced  by  the  wing/ vortex  interaction.  A  general  expression 
for  the  incremental  suction  pressure  is  given  by 


U^b)  (y-y')"  +  (h/b)" 


+  16  —e 


2  -2  (y-y')"  +  h^  /r. 


The  quantity  r  represents  the  strength  of  the  vortex  and  the 

rest  of  the  terms  are  illustrated  in  Figure  3.  The  first 

term  in  Equation  1  represents  the  contribution  of  the  tangential 

velocity  associated  with  the  vortex,  whereas,  the  second  term 

represents  the  effects  of  the  axial  velocity  in  the  core  of  the 

vortex.  This  simplified  expression  was  used  to  obtain  good 

correlation  with  observed  experimental  results  by  selecting 

appropriate  values  of  the  parameters  A  and  h/r  .  Quali- 

c 

tatively,  the  trends  predicted  were  good  but  the  values  of 
A  required  were  somewhat  unrealistic. 

In  Reference  1,  expressions  were  also  derived  for  the 
incremental  lift  2uid  drag  coefficients  due  to  the  effects  of 
line  vortices  interacting  with  a  lifting  surface.  Using  these 
expressions,  comparisons  of  experimental  and  theoretical  results 
were  made  for  configurations  with  a  strake,  a  close  coupled 
Ccmard,  and  a  leading  edge  snag.  These  comparisons  are  shown 
in  Figures  4  through  6  respectively.  As  can  be  seen  from  the 
comparisons,  the  simplified  theory  generally  predicted  the 
overall  change  caused  by  the  vortex  interactions. 

Additional  comparisons  were  also  made  in  Reference  1 
to  assess  the  effectiveness  of  the  different  high  lift  devices. 
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These  comparisons  are  presented  in  Figures  7  through  9.  The 
results  indicated  that  at  low  angles  of  attack,  the  snag  was 
probably  the  most  effective  device,  whereas,  at  high  angles 
of  attack,  the  attached  strake  would  probably  be  the  most 
effective  device  in  improving  the  high  angle  of  attack 
maneuvering  capabilities. 

The  general  conclusions  that  were  drawn  from  the 
preliminary  efforts  summarized  above  from  Reference  1  were: 

a.  The  interaction  of  a  wing  with  free  vortices, 
generated  by  strcOces  or  close  coupled  canards,  gives  rise 
to  beneficial  lift  and  generates  some  drag  increase  at  low 
lift  coefficients  but  lower  drag  coefficients  at  high  lift 
coefficients. 

b.  The  interaction  of  the  wing  with  the  vortex 
generated  within  the  pressure  field  of  the  wing,  such  as  those 
generated  by  leading  edge  vortex  generators  (e.g. ,  snags  and 
strakes),  gives  rise  to  a  relatively  weak  vortex  and  a  bene¬ 
ficial  incremental  lift  through  an  enhanced  suction  field 
much  greater  than  that  produced  by  a  free  vortex  of  the  same 
strength  interacting  with  a  lifting  surface.  In  addition, 

the  incremental  drag  generated  by  the  leading  edge  vortices 
seems  to  be  favorable  throughout  the  range  of  angles  of 
attack  of  interest  to  a  maneuverable  operating  aircraft.  On  the 
basis  of  these  tentative  conclusions  arrived  at  through  the 
use  of  a  simplified  theory,  an  experimental  program  was 
undertaken  to  determine  if,  in  fact,  the  noted  benefits  could 
be  realized. 
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2.2  EXPERIMENTAL  INVESTIGATION 

Based  on  the  initial  theoretical  investigation,  a 
series  of  experimental  investigations  were  undertaken,  to 
demonstrate  and  verify  that  the  beneficial  incremental  lift  can 
be  obtained  from  a  controlled  vortex  generated  by  and  inter-  ^ 

acting  with  a  lifting  surface,  to  determine  the  mechanism  by  ; 

which  the  favorable  performance  gains  were  achieved,  and  to  | 

evaluate  the  effects  of  the  trailed  wake  on  the  effectiveness 
of  a  horizontal  tail  plane.  Studies  were  conducted  on  an  F-4 
type  of  wing  planform  with  and  without  leading  edge  vortex 
control  devices  which  included  a  strake,  a  snag,  a  separator 
plate,  cmd  combinations  thereof.  The  total  forces,  moments, 
and  detailed  pressure  distribution  on  the  lifting  surface  were 
measured  and  flow  visualization  was  carried  out  using  tufts 
and  the  helium  bubble  technique. 

A  systematic  series  of  tests  on  a  basic  wing  planform 
with  and  without  attached  leading  edge  devices  such  as  the 
stredces,  the  snags,  and  separator  plates  were  undertaken  to 
investigate  the  effects  of  attached  vortices  on  the  aerody- 
neunic  characteristics  of  the  lifting  surface.  The  basic  wing 
model  was  a  1/4  scaled  version  of  an  F-4  type  wing  with  a 
NACA  0006  airfoil  section.  The  additional  leading  edge  devices 
were  made  of  flat  beveled  edge  metal  plates  and  could  be  placed 
at  various  locations  on  the  leading  edge  of  the  basic  wing. 

The  line  drawing  presented  in  Figure  10  illustrates  the  manner 
in  which  various  configurations  were  obtained  by  adding  appen¬ 
dages  to  the  basic  lifting  surface.  The  basic  wing  is  shown 
by  solid  lines  and  the  streUce  and  the  snag  are  shown  by  dotted 
lines.  Many  other  configurations  were  formulated  and  tested 
by  this  procedure  but  will  not  be  summarized  in  this  report. 
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Complete  details  on  these  additional  configurations  can  be 
found  in  References  2  and  3. 

The  instrumentation  used  to  record  the  data  during 
these  tests  consisted  of  a  yoke  type  balance  located  beneath 
the  floor  to  measure  the  three  force  emd  the  three  moment 
components.  The  balance  measurements  were  monitored  on-line 
prior  to  offline  computer  processing  into  the  appropriate  non- 
dimensional  components.  The  reference  axes  for  these  measure¬ 
ments  are  shown  in  Figure  11. 

Pressure  measurements  were  made  at  220  tap  locations 
on  the  wing  surface  using  five  48  port  scanivalves,  each  scani- 
valve  monitoring  the  tunnel  total,  static,  and  tunnel  centerline 
pressures.  The  pressure  data  was  directly  punched  on  cards  and 
later  converted  into  pressure  coefficients.  Pressure  and  force 
measurements  were  obtained  at  both  positive  and  negative  angles 
of  attack. 

Flow  visualization  in  the  vicinity  of  the  wing  was  pro¬ 
vided  by  neutrally  buoyeuit  helium  bubbles.  The  bubbles  were 
illuminated  by  a  collimated  beam  of  light  smd  were  observed 
and  photographed.  Flow  visualization  was  also  obtained  using 
tufts  located  at  the  pressure  tap  locations  over  the  surface. 

The  majority  of  the  tests  were  conducted  at  tunnel 
speeds  of  100  and  150  ft/sec  over  an  angle  of  attack  range 
of  zero  to  approximately  thirty  (30)  degrees. 

Since  the  primary  purpose  of  the  wind  tunnel  tests 
was  to  determine  if  the  "vortex  lift"  could  have  a  significant 
effect  on  the  performance  characteristics  of  a  low  aspect  ratio, 
swept  lifting  surface  at  high  angles  of  attack,  the  pressure 
distributions  and  flow  visualization  over  the  wing  were  the 
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primary  tools  used  to  evaluate  the  manner  in  which  the  attached 
vortex  flows  had  altered  the  performance  characteristics  of  the 
lifting  surface  as  measured  by  the  balance  system. 


Detailed  discussion  of  the  results  obtained  during 
these  tests  is  presented  in  References  2  and  3.  Only  the 
major  results  that  were  obtained  will  be  discussed  herein,  with 
the  help  of  representative  data.  A  majority  of  the  character¬ 
istics  of  the  interaction  of  the  wing  and  vortex  can  be  illu¬ 
strated  by  the  discussion  of  the  results  obtained  for  the 
following  three  configurations; 

a.  Basic  wing 

b.  Basic  wing  with  strake 

c.  Basic  wing  with  a  strake  and  a  snag 

The  variations  of  the  lift  coefficient,  the  drag  co¬ 
efficient  and  the  pitching  moment  coefficient  with  angle  of 
attack  for  the  aUDove  configurations  are  presented  in  Figures  12 , 
13  and  14,  respectively.  These  figures  indicate  that  the  per¬ 
formance  characteristics  of  all  three  configurations  are 
approximately  the  same  below  an  angle  of  attack  of  about  13*. 
Beyond  this  angle,  the  basic  wing  shows  a  soft  stall,  with  a 
maximum  lift  coefficient  of  about  1.  The  addition  of  a  strake 
alters  the  lift  curve  considerably  above  13®  angle  of  attack. 

The  lift  coefficient  continues  to  increase  with  angle  of 
attack,  at  approximately  the  same  rate  as  the  one  below  13®, 
up  to  an  angle  of  attack  of  28*  and  then  stalls  abruptly.  At 
an  angle  of  attack  of  28®,  the  strake  configuration  has  a  lift 
coefficient  1. 45  times  the  maximum  value  measured  for  the  basic 
wing.  This  increase  in  lifting  capability  is  due  to  the  strong 
interaction  between  the  strake  and  the  leading  edge  vortices. 

The  addition  of  the. snag  to  the  wing  strake  combination 
had  a  beneficial  effect  in  the  high  angle  of  attack  range  in 
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that  it  increased  the  lift  at  a  given  angle  of  attack  and 
yielded  a  maximum  lift  coefficient  of  about  1.6  times  that 
of  the  basic  wing  and  1.1  times  that  for  the  basic  wing  with 
a  root  strake. 

The  effect  of  the  strake  and  snag  on  the  drag  charac¬ 
teristics  of  the  lifting  surface  are  shown  in  Figure  13.  The 
drag  for  those  configurations  was  greater  than  for  the  basic 
wing,  but  was  not  greater  to  the  degree  that  would  be  expected 
based  on  the  quadratic  variation  of  the  induced  drag  with  lift. 

The  pitching  moment  coefficients  about  quarter  chord 
station  of  the  mean  aerodynamic  chord  (MAC)  are  shown  in  Fig¬ 
ure  14.  As  can  be  seen,  the  vortex  generating  surfaces  have 
a  pronounced  effect  on  the  pitching  moments.  For  the  basic 
wing,  the  nose  down  (negative)  pitching  moment,  increases 
somewhat  linearly  with  the  angle  of  attack  until  the  angle 
of  attack  at  which  leading  edge  separation  occurs  and  the 
attendant  leading  edge  vortex  forms.  At  that  point,  the 
pitching  moment  increases  much  less  rapidly  with  angle  of 
attack.  The  rate  of  increase  of  the  pitching  moment,  how¬ 
ever,  becomes  larger  until  it  reaches  approximately  the  same 
value  as  obtained  before  the  initiation  of  the  leading  edge 
separation.  The  break  in  the  pitching  moment  slope  is  due  to 
the  formation  of  the  leading  edge  vortex.  As  the  angle  of 
attack  is  increased  further,  the  effect  of  the  vortex  diminishes 
due  to  separation  from  the  wing  surface  and  completely  disappears 
at  an  angle  of  attack  of  approximately  22**. 

With  the  strake  and  the  strake  and  snag  added  to  the 
basic  wing,  the  trend  of  the  pitching  moment  with  angle  of 
attack  is  approximately  the  same  below  an  angle  of  attack  of 
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13®  as  that  of  the  basic  wing.  Above  13®,  however,  the  trend 
of  the  pitching  moment  with  increasing  angle  of  attack  is 
opposite  to  that  of  the  basic  wing.  The  reason  for  this 
change  is  that  when  the  leading  edge  vortex  is  formed,  its 
interaction  with  strake  vortex  causes  it  to  be  strengthened 
and  stabilized  and  then  it  continues  to  grow  in  strength  as 
the  angle  of  attack  is  increased.  Since  the  leading  edge  vor¬ 
tex  is  formed  over  the  inner  half  of  the  span,  ahead  of  the 
quarter  chord  point  of  the  MAC,  it  creates  a  strong  nose  up 
moment  and  dominates  the  pitching  moment.  For  example,  at  30® 
angle  of  attack,  the  data  indicates  that  the  pitching  moment 
generated  by  the  leading  edge  vortex  is  1.5  times  greater  than 
that  generated  by  the  entire  lifting  surface.  With  the  snag 
added,  the  leading  edge  vortex  moves  slightly  aft  and  the  ef¬ 
fect  is  comparatively  less  pronounced  than  that  of  the  strake 
wing  configuration. 

Even  though  the  pitching  moment  is  destabilizing  in 
the  post-leading  edge  stall  regime,  its  maximum  magnitude  is 
less  than  one  third  of  that  which  was  measured  for  the  basic 
wing  and,  therefore,  was  thought  to  be  controllable  with  a 
properly  placed  horizontal  tail  surface.  The  proper  placement 
of  the  horizontal  tail,  however,  may  be  significeintly  different 
for  this  configuration  than  it  would  be  for  the  basic  wing  be¬ 
cause  of  the  flow  field  on  the  horizontal  tail  induced  by  the 
vortex  wake  of  the  wing.  It  is  of  interest  to  note  that  since 
the  vortex  lift  dominates  the  pitching  moment  characteristics, 
it  might  be  possible  to  maintain  a  constauit  pitching  moment  at 
all  angles  of  attack  in  the  post  stall  region  by  proper  wing 
planform  tailoring  to  obtain  an  optimum  leading  edge  vortex 
placement. 

The  rolling  moment  characteristics  were  not  altered  sig¬ 
nificantly  because  of  the  addition  of  the  leading  edge  devices. 
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It  was  found  that  the  center  of  drag  tends  to  move  inboard 
and  remain  stationary  with  angle  of  attack,  which, it  is  be¬ 
lieved,  is  advantageous  in  terms  of  yaw  trim. 

No  significant  changes  were  observed  in  the  lift  and 
drag  characteristics  when  the  relative  angle  between  the  strake 
and  the  basic  wing  was  changed.  However,  the  pitching  moment 
characteristics,  as  shown  in  Figure  15,  indicate  a  possibility 
of  obtaining  a  high  degree  of  pitch  control  without  affecting 
the  lift  and  drag  characteristics. 

The  above  performance  characteristics  and  the  manner 
in  which  the  free  vortices  control  the  flow  can  be  explained 
by  the  pressure  distributions  on  the  surface  and  flow  visual¬ 
ization.  Comparisons  of  pressure  distributions  at  13°  (angle 
at  which  the  flow  is  fully  attached)  and  at  27.7°  (angle  at 
which  maximum  lift  was  obtained)  were  made. 

Figure  16  presents  the  spanwise  pressure  distribution 
along  the  constant  chord  lines  for  the  three  configurations 
of  interest  at  an  angle  of  attack  13.1°.  While  the  total  lift 
generated  by  these  configurations  is  almost  the  same  at  this 
angle  of  attack  (Fig.  12) ,  the  pressure  distributions  are 
significantly  different. 

It  can  be  shown  that  these  differences  in  the  pres¬ 
sure  distributions  are  caused  by  the  vortex  flows  that  are 
generated  by  the  lifting  surface  and  the  attached  vortex 
generators. 

In  analyzing  the  effect  of  the  attached  vortex  flows 
on  the  pressure  distribution  over  the  surface  of  the  wing,  both 
the  induced  and  the  suction  effect  of  the  vortex  ilows  must  be 
considered.  While  the  induced  effect  is  generally  generated 
at  spatial  locations  greater  than  the  core  diameter,  the  vortex 
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suction  effect  occurs  at  spatial  locations  associated  with 
dimensions  on  the  order  of  the  vortex  core  diameter. 

The  spanwise  pressure  distributions  presented  in  Fig¬ 
ure  16  for  the  basic  wing  indicate  that  there  is  the  semblance 
of  a  leading  edge  vortex,  which, after  its  formation  as  a  co¬ 
ordinated  vortex  along  the  5%  chord  line,  trails  outward  and 
aft  and  leaves  the  trailing  edge  near  the  wing  tip.  (The 
path  of  the  centerline  of  the  vortex  can  be  traced  by  the 
peaks  in  the  pressure  distribution).  The  pressure  distribu¬ 
tions  also  indicate  that  a  small  secondary  leading  edge  vortex 
is  formed  at  the  leading  edge  at  about  75%  span  location.  This 
secondary  vortex,  which  is  probably  formed  by  the  induced  effect 
of  the  primary  leading  edge  vortex  as  it  turns  in  the  chordwise 
direction,  is  quickly  separated  from  the  surface  of  the  wing  by 
the  induced  upwash  of  the  primary  vortex  flow  from  the  broad¬ 
ness  of  the  peaks  in  the  pressure  distribution  and  the  rather 
low  value  of  the  suction  pressure  of  the  main  leading  edge 
vortex,  it  is  obvious  that  the  leading  edge  vortex  at  this 
angle  of  attack  is  rather  weak  and  diffuse.  This  conclusion 
is  supported  by  the  external  flow  visualization  pictures  that 
were  taken  using  helium  bubbles  as  only  a  slight,  large  diameter 
swirling  flow,  that  might  be  associated  with  a  leading  edge 
vortex,  was  indicated  in  the  pictures.  The  induced  effect  of 
the  leading  edge  vortex,  although  relatively  weak  at  this  angle 
of  attack,  tended  to  reduce  the  angle  of  attack  over  the  inboard 
sections  of  the  lifting  surface  which  in  turn  reduced  the  chord- 
wise  gradient  of  potential  lift  in  this  region. 

The  pressure  distribution  for  the  strake  configuration 
shows  some  overall  similarities  to  that  for  the  basic  wing, 
although,  there  are  differences  in  detail.  The  effect  of  the 
strake  vortex  is  clearly  evident  near  the  root.  The  path  of 
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the  strake  vortex  is  basically  in  the  chordwise  direction. 

Since  the  magnitude  of  the  pressure  peaks  associated  with 
the  streJce  vortex  decreases  as  the  vortex  moves  downstream, 
it  can  be  concluded  that  the  strake  vortex  moves  away  from 
the  wing  as  it  moves  downstream.  Both  of  these  conclusions 
seem  to  be  valid  from  photographs  in  Figure  17.  These  photo¬ 
graphs  are  obtained  from  the  helium  bubble  flow  visualization 
and  indicate  that  the  strake  vortex  goes  directly  across  and 
away  from  the  wing  surface  as  it  moves  downstream.  The  in¬ 
duced  effect  of  the  strake  vortex,  though  relatively  weak, 
creates  spanwise  loading  gradients  near  the  leading  edge  which 
gives  rise  to  the  leading  edge  vortex.  The  induced  effect  of 
the  strengthened  leading  edge  vortex  tends  to  reduce  the  poten¬ 
tial  lift  of  the  inboard  section  of  the  lifting  surface  as  in 
the  case  of  the  basic  wing. 

In  addition,  the  upwash  induced  by  the  leading  edge 
vortex  over  the  outboard  section  of  the  wing  seems  to  separate 
the  flow  as  the  surface  pressures  are  less  than  they  were  for 
the  basic  wing  in  this  region.  Due  to  the  induced  and  suction 
effect  of  the  strake  vortex  and  the  stabilized  flow  over  the 
inboard  sections  of  the  lifting  surface,  the  leading  edge  vor¬ 
tex  not  only  becomes  slightly  stronger  but  much  more  concentrated 
and  thus  the  suction  peak  it  creates  over  the  outboard  sections 
of  the  lifting  surface,  while  narrower,  is  larger  in  magnitude. 

It  can  also  be  seen  by  comparing  the  pressure  distributions  of 
the  wing  with  and  without  the  strake,  that  once  the  leading  edge 
vortex  of  the  wing  strake  configuration  becomes  a  "semi-free" 
vortex  it  is  turned  in  the  free  stream  direction  much  more  rapidly 
than  it  was  for  the  basic  wing  configuration  thus  creating  a 
rather  effective  aerodynamic  wing  fence  across  the  chord  at 
approximately  the  80%  span  location. 
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The  surface  flow  for  this  configuration  as  depicted 
by  tufts  is  shown  in  Figure  18.  The  centerline  of  the  paths 
of  both  the  strake  and  leading  edge  vortices  have  been  in¬ 
dicated  on  those  photographs.  As  indicated  by  the  tufts  the 
flow  in  between  the  streeunwise  positions  of  the  two  vortices 
is  in  the  clockwise  direction  indicating  that  potential  flow 
and  lift  are  being  developed  in  this  area.  While  there  is 
some  effect  of  vortex  induced  flow  inboard  of  the  strake  vortex 
location,  the  flow  is  basically  potential  in  nature  over  that 
region  as  well.  Outboard  of  the  leading  edge  vortex  as  it 
crosses  the  chord,  however,  it  is  noted  that  the  flow  is  basi¬ 
cally  in  the  spanwise  direction.  A  closer  look  at  the  tufts 
along  the  leading  edge  outboard  of  the  leading  edge  vortex 
indicated  that  this  area  is  stalled  and  the  flow  is  separated. 

It  is  believed  that  the  large  induced  effect  of  the  leading 
edge  vortex,  which  creates  an  upwash  in  this  region,  is  respon¬ 
sible  for  the  flow  separation.  Once  the  flow  over  the  surface 
of  the  wing  has  been  separated,  the  induced  velocity  of  the 
leading  edge  vortex  creates  the  spanwise  flow  in  the  low  velocity 
region  of  the  separated  flow.  While  the  pressure  distributions 
presented  for  the  wing  with  both  the  strake  and  outboard  snag 
added  is  considerably  different  than  those  of  the  other  con¬ 
figurations,  the  same  total  lift  is  obtained.  The  pressure 
peaks  associated  with  the  strake,  leading  edge,  and  snag  vortex 
are  clearly  evident  in  the  distributions  that  are  presented  in 
Figure  16.  As  can  be  seen,  the  characteristics  of  the  strake 
vortex  remains  about  the  same  as  they  were  for  the  previous  ’ 
configuration  and  the  leading  edge  vortex  has  been  weakened 
and  moved  inboard  due  to  its  interaction  with  the  rather  strong 
snag  vortex.  Since  the  straJce  vortex  is  rather  far  removed 
from  either  the  leading  edge  or  snag  vortex^  there  is  little 


evidence  of  interacting  effects  of  the  strake  vortex  with  the 
outboard  vortices.  In  fact,  the  pressure  distributions  for 
this  configuration,  except  for  the  independent  effects  of  the 
strake  vortex  are  almost  identical  to  those  of  the  basic  wing 
with  just  the  snag  added.  On  the  basis  of  the  pressure  dis¬ 
tributions  obtained  for  the  various  configurations  at  a  rela¬ 
tively  low  angle  of  attack,  it  is  concluded  that  while  the 
influence  of  the  rather  weak  vortex  flows  did  not  alter  the 
total  lift  that  was  generated,  they  did  alter  significantly 
the  loading  distributions. 

At  angles  of  attack  much  greater  than  the  angle  at 
which  the  basic  wing  stalls,  the  influence  of  the  vortex  flows 
on  the  total  lift  that  is  generated  is  much  more  significant 
as  is  shown  by  the  pressure  distributions  presented  in  Fig¬ 
ure  19.  The  pressure  distributions  shown  for  the  basic  wing 
indicate  that  except  for  a  weak  pressure  peak  due  to  the 
leading  edge  vortex  near  the  wing  root,  the  entire  lifting 
surface  is  completely  stalled  as  indicated  by  pressure  co¬ 
efficients  of  near  unity.  In  contrast,  the  pressure  distri¬ 
butions  obtained  for  the  wing  with  the  strake  added  indicate 
that  large  pressure  peaks  are  obtained.  It  can  also  be  seen 
that  not  only  the  peak  pressure  associated  with  the  leading 
edge  vortex  increased  but  the  position  of  the  vortex  has 
moved  inward.  It  is  believed  that  both  of  these  effects  are 
due  to  the  interaction  of  the  two  vortices  and  the  strengthen¬ 
ing  of  the  leading  edge  vortex  by  the  loading  gradients  gen¬ 
erated  by  the  strake  vortex. 

As  indicated  by  the  pressure  distribution,  the  strake 
vortex  has  an  entirely  different  path  over  the  surface  of  the 
wing  at  27.7  degrees  than  it  did  at  13.1  degrees  angle  of 
attack.  This  difference  in  the  paths  of  the  strake  vertex  can 
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be  seen  by  comparing  the  visualized  strake  vortices  in  Fig¬ 
ure  20  with  those  in  Figure  17.  While  the  initial  rate  of 
movement  away  from  the  surface  of  the  strake  vortex  is  greater 
at  27.7  degrees  than  it  was  at  13.1  degrees,  the  strong  in¬ 
duced  effect  of  the  leading  vortex  moves  the  strake  vortex 
back  towards  the  surface  of  the  wing  as  well  as  outward  along 
the  75%  chordline.  As  shown  in  Figure  21,  the  leading  edge 
vortex  moves  out  along  the  leading  edge  as  well  as  away  from 
the  lifting  surface  due  to  the  induced  effect  of  the  strake 
vortex  until  it  bursts.  After  bursting  the  leading  edge  vor¬ 
tex  becomes  turbulent  and  turns  in  the  streamwise  direction. 
Although  not  shown  in  these  pictures  the  strake  and  burst 
leading  edge  vortex  could  be  seen  to  mix  near  the  trailing 
edge  of  the  lifting  surface  near  the  70%  span  location  as 
they  both  moved  in  the  downstream  direction. 

The  tuft  picture  visualizing  the  surface  flow  for  this 
configuration  at  an  angle  of  attack  at  27.7  degrees  is  shown  in 
Figure  22.  The  small  areas  of  potential  flow  generated  by  the 
induced  effect  of  the  leading  edge  and  strake  vortex  can  be 
seen  in  the  picture.  It  can  also  be  seen  that  the  strong  in¬ 
duced  effect  of  the  leading  edge  vortex  as  it  moves  over  the 
chord  of  the  wing  has  completely  stalled  the  tip  region  and 
generates  a  strong  spanwise  flow  in  the  separated  flow  region 
over  the  mid  span  portion  of  the  lifting  surface  between  the 
two  vortices. 

The  pressure  distributions  presented  for  the  wing  with 
the  strake  and  snag  added  indicate  that  at  this  high  angle  of 
attack  the  primary  effects  of  the  snag  vortex  have  disappeared 
but  the  effects  of  the  strake  and  leading  edge  vortices  have 
remained  about  the  same  as  they  were  for  the  wing-strake 
configuration. 
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It  is  believed  that  the  pressure  distribution  pre¬ 
sented  for  the  various  configurations  at  27.7  degrees  clearly 
indicates  the  strong  influence  of  the  vortex  flows  on  the  lift 
that  is  generated  by  the  wing.  The  majority  of  additional 
lift  obtained  from  the  vortices  is  due  to  the  low  pressures 
within  the  core  of  the  vortices  and  not  due  to  their  induced 
effects  on  the  lifting  surface. 

The  above  discussion  essentially  states  that  a  single 
vortex  formed  within  the  pressure  field  (as  for  the  basic 
wing  alone)  is  rather  unstable  and  ineffective  in  the  post 
stall  region.  In  the  presence  of  another  vortex  (such  as  a 
strake  vortex  created  by  the  root  strake)  there  is  a  mutual 
strengthening  and  stabilizing  that  leads  to  large  increments 
of  lift  well  into  the  post  stall  region. 

It  was  thought  that  the  addition  of  an  outboard  strake 
could  provide  stabilization  of  the  flow  on  the  outboard  sec¬ 
tion  of  the  wing.  A  typical  multi-strake  configuration  is 
shown  in  Figure  23. 

When  compared  with  the  performance  obtained  for  a 
single  strake  configuration,  the  performance  of  such  a  multi- 
strake  configuration  indicated  little  improvement  because  the 
two  free  edges  of  the  outboard  strake  generated  two  vortices 
of  opposite  strength  which  essentially  nullified  the  effects 
of  each  other.  The  induced  effects  of  the  outboard  strake 
vortices  tended  to  weaken  the  main  leading  edge  vortex  which 
is  inboard  of  that  strake.  The  addition  of  the  outboard 
strake,  however,  did  alter  the  pressure  distribution  though 
the  integrated  effect  on  the  performance  is  small. 

Since  the  beneficial  effect  of  the  vortices  were 
realized  at  high  angles  of  attack  and  could  be  traced  to  the 


presence  of  the  leading  edge  vortex,  it  was  thought  that 
similar  benefits  could  be  achieved  at  lower  angles  of  attack 
by  generating  the  leading  edge  vortex  by  premature  separation 
using  a  sharp  edged  separator  plate.  A  configuration  using 
a  strake  and  a  full  separator  plate  is  shown  in  Figure  24. 

The  measured  performance  indicated  that  at  lower  angles  of 
attack,  add' ;.,ional  lift  was  obtained  for  the  configuration 
shown  in  Figure  24,  as  compared  to  the  configuration  with 
only  one  strake.  The  additional  lift,  however,  was  accom¬ 
panied  by  a  significant  increase  in  drag.  The  reason  for 
this  was  that  the  leading  edge  vortex  forms  much  further  in¬ 
board  and  affects  much  larger  areas  of  the  wing  and,  thus,  a 
large  portion  of  the  wing  outboard  was  stalled  due  to  the 
induced  effects  of  the  vortices. 

A  final  configuration  was  tested  that  tended  to  elim¬ 
inate  the  disadvantages  of  the  last  two  configurations.  This 
configuration  is  shown  in  Figure  25.  The  basic  feature  of 
this  configuration  is  that  it  does  not  have  the  additional 
inboard  free  edge  that  characterized  the  single  outboard 
strake.  It  did  have  a  swept  edge  at  the  outboard  end  of 
the  partial  separator  plate  that  provided  the  free  edge  that 
is  desirable  in  order  to  form  a  vortex  at  that  location.  The 
strake  vortex  and  the  leading  edge  vortex  formed  over  the 
partial  separator  plate  combined  together  very  quickly  and 
traversed  the  wing  as  a  single  vortex  and  an  additional  free 
vortex  was  formed  over  the  outboard  region.  With  this  con¬ 
figuration,  larger  lift  could  be  achieved  at  higher  angles  of 
attack  but  with  an  attendant  increase  in  drag. 

The  experimental  research  conducted  so  far  had  been 
concerned  primarily  with  a  1/4-scale  model  similar  to  the 
F4E  wing  planform  and  with  retrofit  leading  edge  vortex  control 
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devices  attached  to  it.  While  the  test  results  on  the  configu¬ 
ration  have  shown  that  positive  performance  gains  can  be  obtained 
by  the  interaction  of  the  wing  with  the  vortices  generated  by 
leading  edge  surfaces,  it  was  found,  however,  that  the  pitching 
moment  tended  to  become  unstable  (nose  up)  at  higher  angles  of 
attack.  This  is  particularly  true  for  configurations  that  uti¬ 
lized  leading  edge  vortex  generators.  This  unstable  pitching 
moment  characteristic  is  realized  as  the  large  incremental  lift 
forward  of  the  MAC  obtained  on  the  inboard  sections  of  the  swept 
wing  is  not  balanced  by  the  lift  generated  by  the  outboard  sec¬ 
tions  of  the  wing  which  are  aft  of  the  pitching  moment  axis  and 
only  carry  small  amounts  of  vortex  lift.  It  was  believed,  how¬ 
ever,  tha€  the  magnitude  of  the  overall  pitching  moment  was  well 
within  the  control  capability  of  a  conventional  tail  plane  oper¬ 
ating  in  undisturbed  air.  Since  the  possible  interaction  of  the 
vortex  flow  generated  by  the  wing  with  the  tail  plane  could 
significantly  reduce  the  effectiveness  of  the  horizontal  tail, 
it  was  concluded  that  the  effectiveness  of  the  horizontal  tail 
operating  in  the  wake  of  the  wing  should  be  investigated.  Be¬ 
cause  of  the  large  effect  that  the  vortices  had  on  the  perform- 
€mce  characteristics  of  the  primary  lifting  surface,  their 
effect  on  the  aerodynamic  characteristics  of  the  tail  plane  was 
investigated  experimentally.  A  detailed  discussion  of  the 
results  obtained  in  these  tests  is  presented  in  Reference  4. 

A  brief  description  of  the  test  results  and  a  summary 
of  the  conclusions  is  presented  below. 

The  model  for  this  experimental  investigation  consisted 
of  two  lifting  surfaces,  a  basic  wing  panel  with  suitable  lead¬ 
ing  edge  attachments  and  a  tail  plane.  Both  the  lifting  sur¬ 
faces  were  semi-span  models  similar  in  geometry  to  and  approxi¬ 
mately  1/4-scale  versions  of  the  F-4  phantom  wing  planform  and 
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its  horizontal  tail.  A  picture  of  the  model  in  the  wind  tun¬ 
nel  is  shown  in  Figure  26. 


The  instrumentation  used  in  these  tests  was  essentially 
the  same  as  used  in  earlier  tests  except  that  additional  pressure 
taps  were  used  on  the  tail  surface  to  monitor  the  pressure. 

The  model  incorporated  provisions  for  mounting  the  tail 
plane  at  various  lateral  and  longitudinal  distances  from  the 
wing  panel  as  well  as  means  of  remotely  controlling  the  angle 
of  attack  of  the  tail  surface  relative  to  the  wing  about  a 
vertical  axis  passing  through  its  mean  aerodynamic  center. 

Each  configuration  of  the  wing-tail  system  tested  dif¬ 
fered  from  the  other  in  two  parameters;  (1)  The  longitudinal 
separation  between  the  wing  MAC  and  tail  MAC  (L) ;  (2)  The 

horizontal  separation  (H)  between  the  wing  mean  aerodynamic 
chord  and  that  of  the  tail.  The  distances,  L  and  H,  were  set 
using  the  turnt2UDle  mechanism  on  the  floor  of  the  test  section. 
The  model  arrangement  on  the  turntable  for  a  typical  configu¬ 
ration  is  shown  in  Figure  27.  In  addition,  the  wing  itself 
had  two  basic  configurations;  (1)  wing  plus  conventional  snag, 
and  (2)  the  wing  plus  snag  with  the  root  strake  added. 

Each  configuration  was  represented  by  a  five  digit 
number  followed  by  a  letter  A  or  B.  The  first  digit  can  be 
either  1  or  2  where  1  refers  to  the  basic  wing  with  an  out¬ 
board  snag  and  2  refers  to  the  basic  wing  plus  snag  and  a  root 
strake.  The  second  and  third  digits  together  give  the  longi¬ 
tudinal  separation  (L)  in  inches.  The  fourth  and  fifth  digits 
together  give  the  offset  (H)  in  tenths  of  root  chord  of  the  tail. 
The  letter  A  denotes  the  wing  alone  configuration  and  the  letter 
B  denotes  the  wing-tail  configuration.  A  typical  representation 
could  be  1-75-06-B,  which  means  a  wing-tail  configuration  with 
an  outboard  snag,  L  »  75  inches,  and  H  »  .06  C^  . 
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Because  of  the  lateral  and  longitudinal  offset  of  the 
horizontal  tail  (Fig.  27),  the  wing  would  be  in  different  loca¬ 
tions  in  the  wind  tunnel  as  compared  to  the  earlier  tests.  The 
basic  aerodynamic  characteristics  of  the  various  wing  configu¬ 
rations  were  measured  prior  to  conducting  the  tests  of  the  wing- 
tail  configurations. 

On  the  basis  of  these  preliminary  tests,  it  was  con¬ 
cluded  that  while  there  were  some  effects  on  the  aerodynamic 
characteristics  of  the  lifting  surfaces  due  to  its  location  in 
the  tunnel,  they  were  not  of  a  magnitude  to  sufficiently  affect 
the  results  of  the  test  program. 

The  Figures  28  and  29  present  the  pitching  moment  vari¬ 
ation  with  the  wing  angle  of  attack  for  the  models  without  the 
tail.  These  results  serve  as  a  reference  for  comparison  with 
the  results  obtained  for  the  models  with  the  tail. 

To  determine  the  effect  of  the  wing  wake  on  the  tail 
effectiveness,  tests  were  conducted  for  the  various  configu¬ 
rations  to  determine  for  various  wing  angles  of  attack,  the 
pitch  euigle  of  the  tail  and  trim,  at  which  the  pitching  moment 
of  the  wing-tail  configuration  went  to  zero.  In  addition,  the 
effectiveness  of  the  tail  surface  in  generating  pitching  moment 
about  the  wing  MAC  was  also  determined  by  varying  the  tail  angle 
of  attack  about  its  trim  value. 

Figure  30  presents  the  results  of  the  tests  conducted 
with  configuration  1-75-06B  which  is  the  basic  wing  configu¬ 
ration  (without  strake) . 

The  pitching  moment  generated  about  the  wing  MAC  by 
the  tail  surface,  as  its  angle  of  attack  was  altered,  was 
surprisingly  linear  regardless  of  the  angle  of  attack  of  the 


main  wing.  The  slope  of  the  pitching  moment  curve  with  tail 
angle  of  attack,  dCj^/da,j,,  provides  a  measure  of  the  effective¬ 
ness  of  the  tail  and  the  results  presented  in  Figure  30  indi¬ 
cate  that  while  there  is  some  loss  in  tail  effectiveness  as 
the  wing  angle  of  attack  is  increased  from  15  to  30  degrees, 
it  is  only  of  the  order  of  25%,  which  is  much  less  than  might 
be  expected. 

The  variation  of  the  tail  angle  of  attack  for  trim 
(relative  to  the  wing  mean  chord)  with  increasing  wing  angle 
of  attack,  shown  in  Figure  30,  indicates  that  above  a  wing 
angle  of  attack  of  20  degrees  very  large  changes  in  the  tail 
geometric  angle  of  attack  are  required  to  maintain  aircraft 
trim.  This  is  explained  below. 

Since  at  lower  angles  of  attack  of  the  wing,  the  tail 
drag  moment  and  the  wing  pitching  moment  tend  to  oppose  each 
other,  the  lift  developed  by  the  tail  would  tend  to  be  small. 
The  pressure  distributions  presented  in  Figure  31  for  the  tail 
when  the  configuration  is  in  trim  for  a  wing  angle  of  attack 
of  8.7  degrees,  indicate  that  indeed  the  tail  lift  force  is 
small. 

As  the  angle  of  attack  of  the  wing  is  increased,  the 
MAC  of  the  tail  surface  will,  with  respect  to  the  free  stream 
direction,  be  located  such  that  its  drag  force  will  create  a 
nose  down  instead  of  a  nose  up  pitching  moment  about  the  wing 
MAC.  For  configuration  1-75-06B  this  shift  occurs  at  a  wing 
angle  of  attack  of  11.3  degrees.  For  wing  angles  of  attack  of 
between  12  and  20  degrees,  the  pitching  moment  of  the  wing  does 
not  increase  very  much  due  to  the  formation  of  the  leading 
edge  vortex  in  this  angle  of  attack  range  (Fig.  28) .  It  would 
be  anticipated  therefore  that  the  tail  angle  of  attack  to 
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maintain  pitching  moment  trim  in  this  angle  of  attack  range 
would  not  have  to  be  changed  significantly. 

As  the  wing  angle  of  attack  is  increased  above  20 
degrees,  the  wing  nose  down  pitching  moment  increases  and 
the  pitching  moment  of  the  tail  drag  force  increases  the 
pitching  moment  about  the  wing  MAC  further.  Thus,  as  the 
angle  of  attack  of  the  wing  increases,  the  drag  moment  arm 
of  the  tail  increases  negatively  and  the  moments  generated 
about  the  wing  MAC  by  the  tail  lift  and  drag  force  tend  to 
offset  each  other.  Since  the  lift  force  increases  more  rapidly 
than  the  drag  force  with  increasing  angle  of  attack,  trim  bal¬ 
ances  can  be  obtained  but  the  angle  of  attack  of  the  tail  sur¬ 
face  must  undergo  a  larger  change  to  maintain  trim  than  when 
the  wing  angle  of  attack  is  lower. 

Comparison  of  the  tail  pressure  distribution  presented 
in  Figure  32  with  those  presented  in  Figure  31  shows  that  the 
tail  surface  is  indeed,  developing  more  lift  at  the  higher 
wing  auigles  of  attack. 

When  the  strake  was  added  to  the  leading  edge,  config¬ 
uration  2-75-06B,  the  pitching  moment  characteristics  were 
different  than  the  basic  wing  (Fig.  28)  but  the  average  tail 
effectiveness  as  a  function  of  wing  angle  of  attack  did  not 
change  to  any  great  degree  (Fig.  33) .  It  is  noted,  however, 
that  the  data  presented  in  Figure  33  indicates  that  there  is  a 
noticeable  variation  in  the  tail  effectiveness  as  a  function 
of  angle  of  attack.  This  variation  seems  to  correlate  with 
the  radical  changes  in  the  tail  geometric  trim  angle  with 

respect  to  the  chordplane  of  the  main  wing.  Inspection  of 
the  pressure  distribution  for  the  wing  strake  -  tail  configu¬ 
ration  shows  why  there  is  a  sufficient  difference  in  the 


38 


behavior  of  the  tail  trim  angle  variation  between  the  two 
configurations.  Figure  34  presents  the  pressure  distribution 
on  the  tail  surface  when  the  wing  tail  configuration  is  in 
trim  at  13.1  degrees.  As  can  be  noted  from  the  pressure  dis¬ 
tribution,  there  is  an  obvious  vortex  interaction  in  the  tail 
surface  at  approximately  the  70%  span.  It  is  believed  that 
this  reversal  of  the  pressure  distribution  in  the  tip  region 
is  caused  by  the  induced  effect  of  the  strake  vortex  crossing 
the  tail  around  65%  of  its  semispan. 

At  a  wing  angle  of  attack  of  21.7  degrees,  no  evidence 
of  induced  effects  from  a  concentrated  vortex  could  be  found 
in  the  pressure  distributions  over  the  tail.  This  can  be 
explained  as  follows.  The  strake  vortex,  under  the  influence 
of  the  strong  leading  edge  vortex  generated  by  the  wing,  moves 
farther  outboard  and  is  located  outboard  of  the  tail  tip. 
Therefore,  the  reversal  of  the  pressure  distribution  near  the 
tip,  as  was  found  at  lower  angles  of  attack,  was  not  present. 
This  is  also  true  for  wing  angle  of  attack  of  30.3  degrees 
where  the  stradce  vortex  is  swept  further  outboard  into  the 
leading  edge  vortex  due  to  its  strong  induced  effect.  As  the 
wing  angle  of  attack  is  increased  above  21.7  degrees,  a  tail 
trim  angle  of  attack  trend  similar  to  that  obtained  with  the 
basic  wing  is  obtained  except  that  the  maximum  change  in  trim 
angle  is  only  5  degrees  instead  of  16  degrees  as  it  was  for 
configuration  1-75-06B.  This  much  smaller  tail  trim  angle 
variation  is  attributed  to  the  fact,  that  for  the  wing  strake 
configuration  at  high  angles  of  attack,  the  pitching  moment 
developed  by  the  tail  drag  and  lift  are  working  together  to 
counter  the  wing  pitching  moment  instead  of  working  against 
each  other  as  they  were  for  the  basic  wing  tail  configuration. 
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The  tail  effectiveness  data  obtained  for  configuration 
1-75-12B  is  presented  in  Figure  35.  These  results  when  com¬ 
pared  to  the  data  presented  in  Figure  30  show  the  effect  of 
moving  the  tail  chordplane  higher  up  from  the  wing  chordplane. 

It  can  be  seen  that  the  results  are  very  similar  for  the  two 
configurations  except  that  for  configuration  1-75-12B,  there 
is  a  sharp  drop  in  the  tail  effectiveness  as  the  wing  angle 
is  increased  from  25.6  degrees  to  29.6  degrees.  It  is  also 
noted  that  the  accompanying  change  in  the  tail  trim  angle  is 
not  as  large  as  it  was  for  configuration  1-75-06B.  It  is 
believed  that  the  reason  the  tail  trim  angle  does  not  vary 
as  much  as  it  did  for  configuration  1-75-06B  in  this  angle 
of  attack  range,  is  that  the  angle  of  attack  at  which  the 
moment  generated  by  the  tail  drag  force  about  the  wing  MAC 
opposes  the  moment  due  to  tail  lift  is  approximately  23  degrees 
instead  of  approximately  11  degrees  as  it  was  for  configuration 
1-75-06B.  Therefore,  the  tail  lift  force  is  much  more  effective 
for  configuration  1-75-12B  to  obtain  trim  than  it  was  for  con¬ 
figuration  1-75-06B. 

It  is  reasoned  that  the  sharp  drop  in  the  tail  effec¬ 
tiveness  as  the  wing  angle  is  increased  from  26  to  approximately 
30  degrees  is  due  to  the  induced  effects  of  the  diffused  leading 
edge  vortex  that  separated  from  the  wing. 

Comparison  of  the  tail  effectiveness  data  presented  in 
Figure  36  for  configuration  2-75-12B  with  that  presented  in 
Figure  33  shows  the  effect  of  doubling  the  distance  of  the 
tail  surface  from  the  wing  chordplane  for  the  wing-strake 
configuration.  It  can  be  seen  that  while  the  tail  effective¬ 
ness  has  a  slightly  different  trend,  the  trend  of  the  tail 
geometric  trim  angle  with  wing  angle  of  attack  is  significantly 
different.  It  is  believed  that  this  significantly  different 
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trend  is  due  to  the  different  orientation  of  the  wing  wake 
with  respect  to  the  tail  plane.  At  21.5  degrees  wing  angle 
of  attack,  the  small  positive  pitching  moment  generated  by 
the  wing  (Fig.  29)  is  countered  by  the  lift  moment  generated 
by  the  tail  surface  and  therefore  the  lift  required  to  be 
generated  by  the  tail  surface  would  be  small.  As  the  wing 
angle  of  attack  is  increased,  the  sudden  jump  in  the  positive 
pitching  moment  generated  by  the  separation  of  the  leading 
edge  vortex  makes  a  demand  for  a  much  greater  lift  from  the 
tail  surface. 

The  effect  of  moving  the  tail  surface  closer  to  the 
wing  trailing  edge  on  the  tail  effectiveness  for  the  wing  and 
wing  strake  configurations  is  shown  in  Figures  37  and  38, 
respectively.  When  this  data  is  compared  to  that  obtained 
for  configurations  1-75-06B,  it  can  be  seen  that  the  results 
and  the  trend  of  the  results  with  the  wing  angle  of  attack 
are  very  similar. 

Figures  39  and  40  present  the  effect  of  the  tail  on 
the  lift  of  the  wing  strake  and  the  basic  wing  configurations. 

It  can  be  seen  from  the  data  presented  that  the  effect  of  add¬ 
ing  the  tail  did  not  significantly  effect  the  total  system 
lift  of  the  wing  strake  configuration.  However,  since  the 
tail  lift  required  to  obtain  pitching  moment  trim  was  in  a 
direction  opposite  to  the  wing  lift  for  the  basic  wing  con¬ 
figuration,  a  significant  reduction  in  the  configuration  lift 
was  obtained  when  the  tail  was  added. 

On  the  basis  of  the  experimental  research,  it  is  gener¬ 
ally  concluded  that  the  test  data  demonstrated  that  an  aircraft 
configuration  can  be  stabilized  with  a  nominal  size  tail  surface 
operating  in  the  vortex  wake  of  the  main  lifting  surface. 
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In  addition  to  the  above  general  conclusion,  the  fol¬ 
lowing  specific  conclusions  were  also  formulated. 

(1)  The  trim  angle  of  the  tail  surface  relative  to  the 
chord  plane  of  wing  varied  differently  with  wing  angle  of 
attack  for  the  basic  wing  and  the  basic  wing  with  the  root 
leading  edge  strake. 

(2)  When  the  concentrated  strake  vortex  from  the  main 
wing  intersected  the  tail  plane,  large  angle  changes  of  the 
tail  were  required  to  maintain  system  trim. 

(3)  Even  though  the  characteristics  of  the  wing  wake  were 

changing  significantly  as  its  angle  of  attack  was  changed  in 
the  range  of  0  ^  30  degrees,  the  tail  effectiveness, 

(dCjjj/da^)  ,  remained  relatively  constant. 

(4)  Even  though  the  tail  surface  could  have  a  geometric 
angle  with  respect  to  the  free  stream  of  24  degrees,  it  is 
estimated,  on  the  basis  of  the  tail  pressure  distributions, 
that  because  of  the  strong  induced  effect  of  the  wing  wake 
that  the  aerodynamic  angle  of  attack  of  the  tail  surface  never 
exceeded  10  degrees  at  trim. 

(5)  On  the  basis  of  conclusions  (3)  and  (4)  it  is  believed 
that  the  total  dynamic  pressure  at  the  leading  edge  of  the  tail 
surface  remained  relatively  constant  over  the  entire  angle  of 
attack  range  of  the  wing. 

(6)  Depending  on  the  longitudinal  and  vertical  location 
of  the  tail  surface  relative  to  the  wing  MAC,  the  drag  of  the 
tail  surface  can  have  a  significant  effect  on  the  pitching 
moment  about  the  wing  MAC. 

(7)  For  the  basic  wing  configuration  having  a  nose  down 
pitching  moment,  the  tail  lift  required  for  system  trim 
tended  to  reduce  the  lift  from  the  value  developed  by  the 
wing  alone.  This  reduction  became  more  significant  as  the 
longitudinal  distance  between  the  wing  and  tail  MACs  was 
reduced. 
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2.3.  THEORETICAL  ANALYSIS 

After  the  experimental  results  obtained  in  the  second 
year  (Ref.  2)  were  analyzed,  it  was  concluded  that  the  simpli¬ 
fied  theory  formulated  in  Reference  1  was  inadequate  to  describe 
the  effect  of  vortex  flows  on  the  low  aspect  ratio  lifting  sur¬ 
faces.  It  was  also  apparent  that  there  was  more  than  one  free 
vortex  above  the  wing  and  at  higher  angles  of  attack  (21.6°,  for 
example) ,  the  flow  was  partially  stalled.  The  method  developed 
in  Reference  1  was  based  on  a  lifting  line  theory  and  could  not 
account  for  the  complex  flow  phenomena  or  the  interactive  effects 
of  one  vortex  on  another.  Therefore,  development  of  a  theoreti¬ 
cal  prediction  technique  which  could  be  applied  to  reasonably 
complex  planforms,  treat  the  multiple  vortices  and  their  inter¬ 
action  with  the  wing  and  each  other  adequately,  treat  the 
separated  flow,  and  have  the  capability  of  obtaining  a  solution 
from  easily  obtained^le  input  data,  was  undertaken. 

Initially,  a  lifting  surface  theory  was  developed.  The 
wing  was  modelled  as  a  surface  of  discontinuity.  The  overall 
prediction  from  this  analysis  was  found  to  be  good,  however, 
certain  problems  were  encountered  in  the  treatment  of  the  sep¬ 
arated  flow.  The  details  of  the  analysis  were  given  in  Ref¬ 
erence  3.  It  was  apparent  that  the  flow  on  the  upper  surface 
was  sometimes  partially  stalled,  although,  the  flow  on  entire 
lower  surface  always  remained  attached.  A  surface  theory  based 
on  satisfying  the  bovindary  conditions  at  the  mean  surface  of 
the  type  developed  in  Reference  3  could  not  account  for  the 
attached  flow  on  one  side  and  separated  flow  on  the  other 
adequately.  It  was  believed,  therefore,  that  the  total  wetted 
surface  of  the  wing  (including  upper  and  lower  surfaces)  must 
be  included  in  the  analysis  if  adequate  prediction  of  the 
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aerodynamic  characteristics  of  low  aspect  ratio  wings  with 
attached  vortex  flows  was  to  be  obtained.  This  three  dimen¬ 
sional  analysis  of  low  aspect  ratio,  highly  swept,  highly 
tapered  wing  with  partial  separation,  interacting  with  the 
free  vortices  will  be  briefly  described  and  the  program  details 
are  presented  in  the  appendix. 

Basic  features  of  such  a  flow  situation  are  schematic¬ 
ally  illustrated  in  Figure  41.  At  high  angles  of  attack,  the 
flow  separates  near  the  leading  edge  and  rolls  into  a  core  of 
rotating  fluid,  thus  forming  a  free  vortex  called  the  leading 
edge  vortex.  The  flow  is  usually  attached  inboard  of  this 
vortex.  Other  free  vortices  such  as  a  strake  vortex  (due  to 
a  strake  upstream  of  the  wing)  and  a  tip  vortex  (due  to  flow 
around  the  wing  tips)  are  also  present  for  most  configurations 
in  addition  to  the  leading  edge  vortex.  The  effect  of  these 
vortices  on  the  flow  over  the  wing  can  be  consider6U3le ,  espe¬ 
cially  on  the  upper  surface;  depending  on  various  configuration 
parameters  and  the  angle  of  attack.  Generally,  the  vortices 
specified  above  induce  downwash  on  the  inboard  section  of  the 
wing  and  upwash  on  the  outboard  section  which  can  result  in 
a  region  of  separated  or  stalled  flow  over  the  outer  wing 
panels.  The  region  where  separated  flow  may  be  encountered  is 
indicated  in  Figure  41. 

These  free  vortices,  as  was  found  in  References  1,  2 
and  3  can  produce  a  strong  negative  pressure  region  in  the 
upper  surface  of  the  wing.  It  has  been  assumed  in  the  analy¬ 
sis  that  this  vortex  suction  effect  can  be  determined  sepa¬ 
rately  from  the  other  effects  on  the  pressure  at  the  wing 
surface.  The  pressure  at  any  point  on  the  wing  is,  therefore, 
considered  to  be  the  sum  of  the  suction  pressure  due  to  free 
vortices  and  the  pressure  due  to  other  effects  such  as  angle 
of  attack,  thickness,  camber,  separation,  etc. 
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In  the  attached  and  separated  region,  therefore,  fol¬ 
lowing  relations  were  taken  to  be  valid 


total 


'Patt  * 


for  attached  region 


and  (2) 

Cp  “  separated  region 

total  sep 

where,  C-  is  the  total  pressure  coefficient, 

total 


C_  is  the  pressure  coefficient  in  the  attached 
att 

flow  region  due  to  all  the  effects  except 
suction 


ACp  is  the  incremental  suction  pressure  coefficient, 
and, 

Cp  is  the  pressure  coefficient  in  the  separated 

S6P 

region  due  to  all  the  other  effects  except 
suction. 


The  task  of  the  theoretical  analysis  was  to  predict 

C  and  ACp  . 
sep  V 


Once  the  previously  noted  assumption  regarding  the 
vortex  suction  effect  is  made,  the  theoretical  representation 
can  be  divided  in  a  similar  manner.  The  suction  effect, 

ACp  can  be  calculated  at  points  on  the  wing  by  applying  the 

V 

self-similar  solution  of  the  Navier-Stokes  Equations.  The 

potential  pressure,  Cp  ,  is  computed  by  solving  the  potential 

*^att 

flow  over  the  attached  region.  An  empirical  method  is  used  to 

determine  the  separated  region  on  the  wing  surface  and  the 

pressure  there,  Cp  ,  is  calculated  using  a  modified  form  of 

sep 
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Bernoulli's  equations.  The  flow  representations  will  be 
discussed  individually  in  the  following  sections. 


2.3.1.  Representation  of  the  Attached  Flow  Region 


The  flow  over  the  attached  region  is  considered  to 
be  potential.  The  Laplace  equation,  in  terms  of  the  velocity 
potential,  is  solved  with  flow  tangency  condition  at  the  wing 
surface  and  zero  disturbance  condition  at  infinity.  A  solu¬ 
tion  to  this  problem  can  be  written  in  terms  of  a  doublet 
density  distribution  over  the  wing  surface  with  the  doublet 
axes  pointing  in  the  outward  normal  direction  on  the  wing 
surface.  The  doublets  automatically  satisfy  the  Laplace 
equation  and  the  infinity  boundary  condition.  The  strength 
distribution,  D,  is  obtained  by  satisfying  the  tangency  con¬ 
dition  on  the  wing  surface.  The  flow  tangency  condition  at 
the  wing  surface  can  be  written  as  follows : 


J  D(?) 


n(?) 


IRf 


dS 


n(rp) 


0  (3) 


where  r  is  the  position  vector  of  a  point  on  the  surface 
D (r)  is  the  doublet  distribution  density  at  a  point 
on  the  surface 

n (r )  is  the  outward  unit  normal  vector  at  point  r  on 
the  surface 

Fp  is  the  position  vector  of  the  point  on  the  surface 
where  the  boundary  condition  is  applied 
R  »  (r-Fp) 

A 

n (F  )  is  the  normal  unit  vector  in  the  outward  direc- 
P  _ 

tion  at  the  point  given  by  rp 

Uq  is  the  free  stream  velocity  vector,  and 
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qv(rp)  is  the  externally  induced  velocity  at  the  point 

7  on  tne  surface 
P 

The  quantities  used  in  this  equation  are  referred  to  a  body 
fixed  coordinate  system  as  shown  in  Figure  42.  Analytical 
solutions  of  Equation  3  for  the  geometry  of  the  wing  under 
consideration  are  impossible,  therefore  Equation  3  is  solved 
numerically.  The  wing  is  divided  into  a  specified  number  of 
panels  over  which  the  doublet  density  is  assumed  to  be  constant. 
Equation  3  then  can  be  written  as 

NA 

T  a.,  d.  »  w  i  =  1,2  •••  N,  (4) 

^  i  j  j  n  j  A 

]»i  ^ 

where  is  the  number  of  panels  in  the  attached  region  of 

the  wing 

is  the  inward  normal  velocity  induced  by  a  doublet 
panel  j  of  unit  density  at  a  control  point  on 
panel  i  ' 

d^  is  the  value  of  doublet  density  on  panel  j,  and 

w^  is  the  outward  normal  velocity  due  to  vortices 
^  and  free  stream  at  the  control  point  on  panel  i 

Figure  42  indicates  a  typical  semiwing  and  geometric 
panel  arrangement  on  the  wing  surface  (note  that  these  panels 
lie  on  the  wetted  surface  of  the  wing  and  not  on  the  mean 
surface) .  The  centroid  of  a  typical  panel  is  selected  as 
its  control  point. 

The  computation  of  is  carried  out  by  using  the 
following  formula: 
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n(r)  »  R 
1ft’ 


dS 


n  (r^) 


where,  IT  *  (r  -  ,  and 

Sj  refers  to  the  area  of  the 


element 


Computation  of  is  facilitated  by  replacing  the  uniform 

density  doublet  panel  with  its  exact  potential  equivalent, 
a  vortex  ring  of  the  same  strengtn  along  the  boundary  curve 

enclosing  the  doublet  area.  The  velocity  influence  coef¬ 
ficients  are  computed  using  the  Biot-Savart  law  for  each 
segment  of  the  boundary  vortex  ring.  It  should  be  noted 
that  due  to  symmetry,  an  exactly  similar  panel  (box)  exists 
on  the  other  half  of  the  wing  with  same  strength  dj.  The 

therefore,  must  include  the  effect  of  this  box  as  well 


as  the  one  shown  in  Figure  30. 


For  a  wing,  there  exists  a  wake  vortex  sheet  behind 
the  trailing  edge.  The  wake  vortex  sheet  can  also  be  rep¬ 
resented  by  a  distribution  of  uniform  density  doublets.  The 
geometry  of  the  wake  sheet  is  prescribed  in  this  analysis. 
Since  the  wake  sheet  is  iniated  along  the  bisector  of  the 
trailing  edge  angle,  the  trailing  edge  Kutta  condition  is 
satisfied. 


The  strength  of  the  doublet  panels  in  the  wake  can 
be  related  to  the  strength  of  the  doublets  at  the  trailing 
edge  of  the  wing  via  the  condition  that  there  be  no  pressure 
discontinuity  across  the  wake  surface.  For  a  steady  flow, 
this  dynamic  condition  can  be  represented  by  specifying  that 
the  doublet  strength  discontinuity  is  constant  along  the  wake 
in  the  downstream  direction.  Quantitatively,  the  strength  of 
a  wake  panel  can  be  related  to  the  difference  of  doublet 
strengths  at  the  trailing  edge  of  the  wing. 
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t  .e. 


(6) 


where  d  is  the  density  of  doublets,  and  the  subscripts  refer 
to  the  locations  of  the  doublet  panels. 

The  influence  coefficients,  for  the  doublet  boxes 

at  the  trailing  edge  can  be  adjusted  to  include  the  appropriate 
contribution  from  the  wake  element  adjacent  to  such  a  trailing 
edge  box.  A  typical  wake  panel  used  in  the  analysis  is  indi¬ 
cated  in  Figure  42  along  with  its  vortex  equivalent.  This 
wake  panel  also  has  its  image  on  the  other  half  of  the  wing 
and  the  effect  of  the  image  is  included  in  the  computation  of 
the  velocity  influence  coefficient.  Presently,  the  wake  ele¬ 
ments  are  assumed  to  extend  to  infinity  downstreeun. 

The  righthand  side  of  Equation  4 ,  w^^  ,  is  the  outward 

normal  velocity  at  a  control  point  on  box  i,  and  it  includes 
the  component  of  the  free  stream  as  well  as  the  induced  con¬ 
tribution  from  the  free  vortices.  It  is  given  by  Equation  7 
below 


w_ 


Uq  Sin 


(Uo  -  e.)  - 


(w.,  ) 


n 


V  1 


(7) 


The  extent  of  the  wing  surface  over  which  Equation  3 
is  applicable  (the  attached  region)  is  determined  in  an  em¬ 
pirical  manner  (to  be  described  later) .  The  unknown  d^  are 

assigned  to  the  boxes  in  the  attached  region  only  and  Equa¬ 
tion  4  is  applied  at  the  control  points  of  those  boxes  only. 
The  resulting  set  of  equations  has  as  many  unknown  d^  as  the 
number  of  equations  and  can  be  solved  for  the  unknown  d^s. 

Once  these  doublets  strengths  are  known,  the  induced 
velocity  due  to  the  doublets  at  any  point  on  the  surface  can 
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be  calculated.  The  vortex  ring  equivalents  used  in  the  cal¬ 
culation  of  the  normal  velocity  are  used  in  the  calculation 
of  the  total  velocity  also.  However,  in  the  computation  of 
the  tangential  velocity  on  the  wing  surface,  care  must  be 
teiken  to  account  for  the  singularity  in  the  computation  of 
effects  of  the  doublets.  This  is  accomplished  by  calculat¬ 
ing  local  velocity  at  a  control  point  on  the  surface  using 
the  following  formula 

\  ^  ^  ^ 

(q^)  -  Vd./2  +  Z  (a^j  i  +  b. j  j  +  c.j  k)  d.  (8) 

i+j 

where  Vd^  is  the  gradient  of  the  doublet  density  strength 
at  a  control  point  on  the  box  i 
a^j,bij,Cij  are  the  velocity  influence  coefficients  due  to 

doublets  of  unit  density  on  box  j  at  a  control 
point  on  the  box  i  in  x,y,z  coordinate  directions 
respectively 

is  the  velocity  due  to  the  doublets 

This  formula  is  obtained  by  taking  the  limit  of  the  integral 
such  as  the  one  in  Equation  5  when  i  =  j.  The  velocity  com¬ 
ponents  such  as  a^j  etc.  are  obtained  by  using  the  vortex 

equivalents  of  the  doublet  panels. 

The  total  velocity,  q^,  is  obtained  by  adding  the  free 

streaun  velocity  to  (q^j)  .  The  pressure  coefficient  at  a  control 

point  on  the  box  i  is  given  by  the  Bernoulli ' s  equation  as 


(C. 


att 
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u. 


where  is  the  free  stream  speed 


(9) 
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2.3.2.  Pressure  Coefficient  in  the  Separated  Region 

The  determination  of  the  separated  region  is  made 

based  on  a  comparison  of  the  aerodynamic  angle  of  attack 

(the  geometric  angle  of  attack  minus  the  induced  angle  of 

attack)  and  an  empirical  stall  angle,  a  .  Actual  compari- 

s 

son  in  the  computer  program  is  made  between  ,  the  outward 

normal  velocity,  and  a  quantity  U  Sin  a_.  The  box  is  con- 

o  s 

sidered  separated  if 

Wn.  >  Sin  Og  (10) 

outward 

It  is  recognized  that  the  computation  of  the  pressure 
coefficient  in  the  separated  region  is  difficult  because  of 
the  complex  nature  of  the  flow.  An  expression  based  on  the 
Bernoulli  equation  is  used  to  determine  the  pressure  coef¬ 
ficient  in  the  separated  region.  The  expression  for  Cp 
is  given  below 


(Cp  )  a  1 

sep  ^ 


AHi 


(11) 


where  q^  is  total  velocity  at  the  control  point 
of  the  separated  box,  and 

is  the  loss  in  the  total  pressure  due  to  viscosity 

The  loss  in  the  total  pressure  (AH)  is  usually  small  and  hence 
is  neglected  in  the  computation  of  Cp 

sep^ 
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2.3.3. 


Pressure  Coefficient  Due  to  Vortex  Suction 


As  mentioned  before,  the  effects  of  the  free  vortices, 
such  as  the  leading  edge  vortex,  the  strake  vortex,  or  the  tip 
vortex,  are  computed  separately.  The  interaction  of  the  free 
vortex  with  the  wing  is  separated  into  the  two  basic  mechanisms 

(1)  the  influence  of  the  vortex  induced  velocity  field,  and 

(2)  the  suction  effect  of  the  viscous  core  on  the  near  field 
pressure  distribution. 

The  induced  effect  is  included  in  the  potential  solution  via 
the  term  |w^  j  .  The  manner  in  which  the  suction  effect  is 

considered  will  be  outlined  in  the  following  section. 

The  near  field  pressure  effects  of  the  viscous  cores 
of  the  free  vortices  are  obtained  by  using  a  set  of  self¬ 
similar  solutions  of  the  Navier-Stokes  equations.  The  coor¬ 
dinate  system  and  the  typical  models  used  are  shown  in  Fig¬ 
ure  43.  The  details  of  the  self-similar  solutions  and  their 
applications  to  the  free  vortices  in  this  problem  are  given 
in  Reference  4.  A  brief  outline  of  the  analysis  is  presented 
here. 

Basically,  two  viscous  core  models  are  used  for  each 
of  the  free  vortices.  In  the  region  near  the  point  of  gene¬ 
ration  (o  <  s  <  z^)  in  Figure  43,  the  core  radius  of  the 
vortex  is  assumed  to  vary  linearly  with  distance  along  the 
centerline,  s.  (The  parameter,  z^,  is  given  as  an  input). 

In  this  region,  the  axial  velocity  in  the  core  is  much  larger 
than  the  outside  velocity.  The  self-similar  solutions 
developed  by  T.R.  Goodman  in  Reference  1  are  used  to  obtain 
the  incremental  suction  pressures  due  to  the  vortex  in  terms 
of  the  core  radius,  vortex  strength  and  other  local  parameters. 
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The  formulas  obtained  in  Goodman's  model  are  given  below 


where  r^  is  the  core  radius 
c 

a  is  the  turbulent  eddy  viscosity  parameter 

s  is  the  distance  along  the  centerline  of  the  core 

U  ,  are  the  axial  and  tangential  velocities  at  the 
in2ix 

W  center  of  the  core  and  periphery  respectively 
insix 

r  is  the  circulation  strength  around  the  vortex 
is  turbulent  eddy  viscosity 

ACp  is  the  incremental  suction  pressure  coefficient 
^  due  to  free  vortices 

d  is  the  distance  of  a  point  from  the  center  of 
the  core  in  radial  direction 
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Uq  is  the  free  stream  speed,  and 

^u'^w'^r'^1'^2'^3'^4'^5  constants  which  are 

obtained  from  the  numerical  solution. 


The  turbulent  eddy  viscosity  parameter  is  given  by 


a  -  +  .6  for  10’  <  R  <  5  X  10’ 


where  is  the  Reynolds  number  and  is  given  by 


=  Hr 


The  actual  contribution  due  to  different  vortices  are  appro¬ 
priately  added  to  get  the  total  suction  effect. 

The  axial  velocity  in  the  core  of  the  vortex  decreases 
with  distance  s  and  at  a  point  downstream  becomes  roughly 
equal  to  the  free  stream  speed.  (This  point  is  provided 
through  the  parameter  as  an  input) .  Beyond  this  point, 

Batchelor's  vortex  core  model  is  used.  The  radius  in  this 
region  (s  >  z^)  is  assumed  to  be  proportional  to  the  square 

root  of  distance  from  a  virtual  origin.  Self-similar  solutions 
developed  by  Batchelor  are  used  in  this  region.  The  formulas 
obtained  from  this  model  are  given  below 


=  - - -  “  ^ 

^  2tt  (d/r^,)  r^  \ 


-M^d* 

e 
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and 


2 


+  V 


where 

o 


(13) 


r 


2tt 


V^,V^  are  the  tangential  and  axial  velocities  at  a  point 
"d"distance  away  from  the  core  center  in  the  radial 
direction 


are  constants 


The  numerical  values  of  the  constants  used  in  the  above 
equations  are  given  in  the  documentation  in  the  appendix. 

In  general  the  strength  eind  the  geometry  of  the  vortex 
is  not  known  apriori.  The  estimates  of  the  strength  are  made 
from  the  leading  edge  bound  vorticity  strength  obtained  in  the 
totally  potential  solution  (no  vortices,  no  separation).  The 
geometry  of  the  free  vortices  is  obtained  by  making  it  force 
free  under  the  prevailing  velocity  field,  that  is,  the  tangent 
to  the  centerline  of  the  vortex  is  made  parallel  to  the  velocity 
vector  at  that  point  on  the  vortex.  The  velocity  vector  under 
consideration  here  is  given  by 

q  “  ^v  ^a  ^o 

where  q  is  the  total  velocity  vector 

q^  is  the  velocity  induced  by  wing  doublets 
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is  the  velocity  induced  by  the  vortices  including 
the  self-induced  velocity 

is  the  axial  velocity  present  in  the  core  of  the 
vortex 

Uq  is  the  free  stream  velocity. 

It  is  obvious  at  this  point  that  velocity  q  in 
Equation  14  depends  on  the  wing  doublet  solution  and  velocity 
q^  in  Equation  3  depends  on  the  strength  and  the  geometry 

of  the  free  vortices.  The  combined  problem  of  obtaining  the 
doublet  strength  on  the  wing  and  the  geometry  and  the  strength 
of  the  free  vortices  is  a  nonlinear  one.  Therefore,  an  itera¬ 
tive  procedure  is  used,  wherein,  first,  a  geometry  and  strength 
for  the  vortices  is  prescribed  as  an  initial  guess.  The  doublet 
strength  on  the  wing  is  computed  based  on  this  input  geometry. 
The  vortex  geometry  is  then  displaced  to  satisfy  the  force  free 
condition.  If  the  displaced  geometry  is  significantly  different 
from  the  initial  geometry,  the  process  of  obtaining  the  doublet 
solution  emd  new  vortex  geometry  is  repeated.  The  final  solu¬ 
tion  is  obtained  when  there  is  no  significant  difference  between 
the  two  consecutive  geometries  of  the  free  vortices. 

The  pressure  calculations  are  done  only  after  the  itera¬ 
tive  process  discussed  above  is  completed. 

2.3.4  Theoretical  and  Experimental  Correlation 

A  computer  program  based  on  the  analysis  presented 
earlier  was  developed  and  used  to  compute  the  loads  and  pres¬ 
sure  distribution  on  the  Navy  F-4  type  wings.  Several  appli¬ 
cations  of  the  program  to  the  basic  wing  and  basic  wing  and 
strake  configurations  were  presented  in  Reference  4. 
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The  program  used  in  the  above  computations  was  stream¬ 
lined  to  facilitate  its  use  by  outside  users.  Some  of  the 
changes  involved  making  the  coordinate  system  consistent  with 
that  used  by  others  in  this  field.  Certain  input-output  for¬ 
mats  were  changed  to  clarify  the  usage.  Also,  certain  repe¬ 
titious  calculations  were  eliminated  resulting  in  a  somewhat 
more  efficient  execution  of  the  program.  A  documentation  of 
the  program  was  prepared  and  is  included  in  this  report  as  an 
appendix. 

The  validity  of  the  potential  flow  model  used  in  the 
program  was  demonstrated  in  Reference  4.  The  program  was  used 
to  compute  the  pressures  and  the  loads  on  the  wing  at  13.1® 
angle  of  attack  with  the  assumptions  that  there  are  no  free 
vortices  in  the  flow  field  and  the  flow  over  the  entire  wing 
is  fully  attached. 

The  pressure  distribution  along  a  chord  located  at 
approximately  15%  semi-span  is  compared  with  experimental 
values  in  Figure  44.  This  spanwise  station  was  selected 
because  it  was  thought  that  the  influence  of  a  relatively  weak 
leading  edge  vortex  is  negligible  at  this  station.  These 
results  indicate  a  fairly  good  correlation  except  near  the 
trailing  edge  where  some  numerical  problems  could  be  expected 
because  of  the  relative  closeness  of  the  upper  and  lower  sur¬ 
faces.  The  value  of  the  total  integrated  lift  coefficient  cal¬ 
culated  was  .65  which  agreed  well  with  the  experimental  value 
of  .  68 . 

The  program  was  then  used  to  compute  pressures  and 
loads  on  the  Navy  wing-strake  combination  at  higher  angles  of 
attack  of  21.6®  and  27.7®.  The  wing  planform  and  the  geometry 
of  the  free  vortices  is  shown  in  Figure  45  for  the  angle  of 
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attack  of  21.6“.  The  geometry  of  the  free  vortices  was  ob¬ 
tained  from  the  flow  visualization  results  and  pressure  dis¬ 
tributions  obtained  experimentally.  The  vortices  were  made 
to  follow  the  paths  of  the  pressure  peaks  observed  on  the  wing. 

Certain  changes  were  made  in  the  input  to  the  program 
as  compared  to  the  input  used  in  the  runs  made  to  obtain  re¬ 
sults  presented  in  Reference  4.  The  box  structure  was  made 
identical  on  both  the  upper  and  lower  surface  with  96  boxes 
representing  each  surface.  Also,  the  initial  portion  of  the 
leading  edge  vortex  was  modeled  to  include  elements  with  vari¬ 
able  vortex  strengths.  This  change  was  made  to  better  represent 
the  growth  in  the  vortex  strength  in  the  region  where  vortex 
is  being  fed  from  the  wing.  The  vortex  elements  following  this 
variable  strength  region  are  assumed  to  have  a  constant  streng^a. 

The  formation  and  the  growth  of  the  leading  edge  vortex 
takes  place  in  the  leading  edge  region  of  the  wing  and,  there¬ 
fore,  its  strength  is  derived  from  the  doublet  strength  solu¬ 
tion  at  the  control  points  on  the  leading  edge  boxes.  The 
variation  of  the  difference  in  the  doublet  strengths  on  the 
upper  and  lower  surface  at  control  points  on  the  leading  edge 
boxes  is  used  to  obtain  the  strength  of  the  variable  strength 
elements  of  the  vortex.  The  vortex  is  assumed  to  be  fed  from 
the  wing  until  the  point  where  the  above  doublet  strength  dif¬ 
ference  variation  at  the  leading  edge  reaches  a  maximum.  There¬ 
after  the  vortex  is  assumed  to  have  a  constant  strength  equal 
to  the  maximum  of  the  difference  in  doublet  strength  at  the 
leading  edge. 

The  comparison  of  the  predicted  pressure  coefficients 
with  the  experimental  values  is  shown  in  Figures  46a  through 
46f.  These  plots  compare  the  pressure  coefficients  at  several 
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spanwise  locations  on  certain  constant  x/c  lines.  It  is  seen 
that  the  agreement  between  the  calculated  and  experimental 
results  is  fairly  good.  Some  discrepancy  exists  near  the  tip 
region  and  in  the  aft  region  of  the  wing  (for  x/c  =  .6  and  .8). 

Exaunination  of  the  predicted  results  indicates  that 
this  discrepancy  is  due  to  the  pressure  computation  in  the 
separated  boxes  and  in  the  attached  boxes  immediately  next 
to  the  separated  region.  Since  the  doublet  strength  in  the 
separated  boxes  is  forced  to  be  zero,  the  velocities  computed 
at  the  control  points  on  the  attached  boxes  are  unduly  large 
making  the  pressure  coefficients  there  unrealistic.  Also,  the 
actual  doublet  distribution  in  the  attached  region  gets  con¬ 
siderably  altered  when  separated  region  exists  on  the  wing. 
This  distribution  is  particularly  affected  near  the  separated 
region  which  acts  as  a  'hole"  in  the  potential  body  since  no 
no-flow  boundary  conditions  are  satisfied  at  the  boxes  in  the 
separated  region. 

It  is  believed  that  inclusion  of  the  separated  boxes 
in  such  a  manner  is  probably  not  accurate  and  yields  unrealis¬ 
tic  results.  An  alternative  approach  to  this  problem  can  be 
suggested  in  which  the  flow  situation  is  initially  considered 
to  be  entirely  potential  (no  separation)  so  that  no-flow  bound¬ 
ary  condition  is  satisfied  over  the  entire  wing.  The  resulting 
solution  is  then  used  to  yield  the  pressure  distribution  which 
is  used  to  determine  the  separation  region  and  the  pressure  in 
that  region.  Certain  criteria  based  on  the  impressed  potential 
pressure  distribution  and  the  boundary  layer  theory  will  have 
to  be  developed  for  this  purpose. 


The  program  is  also  seen  to  overpredict  the  pressure 
on  the  lower  surface  at  x/c  =0.05.  It  is  thought  that  the 
leading  edge  separation  occurs  on  the  lower  side  of  the  wing 
considering  the  high  angle  of  attack  whereas  the  analysis  con¬ 
siders  the  whole  lower  side  of  the  wing  to  be  attached.  Also,  a 
certain  cunount  of  numerical  inaccuracies  exist  in  this  region 
since  the  chordwise  size  of  the  boxes  in  this  region  is  very 
small. 

In  terms  of  the  suction  peaks  and  their  location,  the 
agreement  seems  to  be  reasonably  good.  The  peaks  are  still 
(as  compared  with  the  results  shown  in  Reference  4)  slightly 
too  narrow  and  too  large,  however,  the  error  can  be  considered 
to  be  small.  This  agreement  can  be  improved  by  further  adjust¬ 
ing  the  vortex  parameters  such  as  the  core  radius  of  the  vortex. 

The  wing  planform  and  the  final  geometry  of  the  vor¬ 
tices  used  in  the  computation  for  the  angle  of  attack  of  27.7® 
is  shown  in  Figure  47.  The  vortex  geometry  in  this  case  is 
also  obtained  from  the  flow  visualization  pictures  and  experi¬ 
mental  pressure  distributions.  The  comparison  of  the  calculated 
and  observed  results  is  shown  in  Figures  48a  through  48d.  It 
can  be  seen  that  the  agreement  between  the  computed  and  observed 
results  is  fairly  good  near  the  forward  portion,  whereas,  the 
agreement  deteriorates  in  the  aft  portions.  The  reasons  for 
these  discrepancies,  as  suggested  earlier,  may  be  due  to  the 
inaccurate  modeling  of  the  separated  flow.  The  separated  flow 
effects  are  probably  stronger  at  this  angle  of  attack  than  at 
21.6®* 

The  integrated  lift  performance  of  the  wing  at  angles 
of  attack  at  13.1®,  21.6®  and  27.7®  is  shown  in  Figure  49.  It 
is  seen  that  the  total  lift  coefficient  in  all  three  cases 
agrees  very  well  with  the  observed  values. 
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It  is  noted  that  these  results  are  based  on  certain 
empirical  parameters  such  as  the  initial  core  radius,  the 
vortex  strengths  of  the  straJce  and  tip  vortices,  and  the 
stall  angle  used  in  the  determination  of  the  separated  region. 
With  a  good  choice  of  these  parameters  the  results  obtained  in 
the  above  cases  appear  to  be  quite  good. 

The  present  program  was  also  used  to  compute  the  pres¬ 
sure  distribution  for  a  delta  wing.  The  delta  wing  configuration 
for  which  the  calculations  were  made  is  shown  in  Figure  50  with 
the  vortex  geometry  superimposed.  The  geometric  data  used  here 
was  obtained  from  Reference  8.  The  delta  wing  under  consider¬ 
ation  is  a  flat  delta  wing  with  aspect  ratio  1.15  and  beveled 
sharp  edges.  The  flow  on  this  wing  separates  at  the  leading 
edge  almost  immediately  as  the  angle  of  attack  is  increased 
from  zero  and  forms  a  vortex  sheet  that  eventually  rolls  up 
into  the  leading  edge  vortex.  An  estimate  of  the  vortex 
geometry  of  the  rolled  up  vortex  core  was  done  based  on 
Smith's  conical  theory.  The  strength  of  the  vortex  based  on 
the  doublet  solution  in  the  leading  edge  region,  however, 
yielded  rather  low  values.  A  probable  reason  for  this  may  be 
that  the  development  of  the  vortex  in  this  case  is  a  much 
slower  process  as  compared  to  that  in  the  Navy  wing  and  should 
in  fact  include  a  roll-up  sheet  which  is  presently  not  included 
in  the  program.  A  somewhat  empirical  approach  was,  therefore, 
taken  to  obtain  results  for  the  delta  wing.  Using  geometry 
obtained  from  Smith's  conical  theory  and  by  adjusting  the 
vortex  strength,  a  good  agreement  was  obtained  between  the 
predicted  total  lift  and  the  experimental  value.  The  computed 
pressures  for  this  case  are  presented  in  Figures  51a  through 
51c.  The  figures  show  the  plots  of  the  pressure  coefficients 
along  the  constant  x  lines.  These  comparisons  indicate  that 
the  program  underpredicts  the  pressure  in  the  root  section  near 
the  apex  (e.g.  x/Cj^  =  .4)  and  it  also  predicts  large  pressure 
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pe^U^s  very  near  the  leading  edge  of  the  wing.  The  under¬ 
prediction  of  pressure  near  the  root  is  probably  because  the 
value  of  the  core  radius  used  in  this  case  is  smaller  than 
the  one  that  should  be  used  since  the  vortex  over  a  delta 
type  wing  is  less  tightly  rolled  compared  to  that  over  a  low 
aspect  ratio  swept  wing.  The  peaks  near  the  leading  edge  are 
the  result  of  the  leading  edge  singularity  which  results  from 
not  satisfying  the  leading  edge  separation  (the  Kutta  condition) 
condition  and  also  from  the  lack  of  the  presence  of  a  roll-up 
vortex  sheet.  It  is  noted  here  that  the  entire  flow  on  the 
delta  wing  was  forced  to  be  attached. 

It  is  noted  that  in  all  the  computations  presented 
above,  a  good  knowledge  of  the  vortex  location  to  be  inputted 
to  the  program  is  necessary.  It  is  very  difficult,  under  the 
framework  of  present  theory,  to  obtain  a  stable  force-free 
geometry  starting  from  an  arbitrary  starting  geometry.  The 
procedure  used  to  obtain  the  force-free  geometry  aligns  the 
geometry  of  the  present  vortex  with  the  calculated  total 
velocity  and  the  calculated  total  velocity  tends  to  pull  the 
vortex  to  the  free  stream  direction.  This  results  from  the 
fact  that  the  induced  velocity  at  the  control  points  on  the 
vortex  is  not  sufficient  to  keep  the  vortex  from  moving  in 
the  free  stre£un  direction.  As  a  matter  of  fact,  the  spanwise 
velocity  at  certain  control  points  appear  to  aggravate  the 
situation  further  instead  of  alleviating  it. 

The  movement  of  the  vortex  in  the  vertical  direction 
is  somewhat  more  orderly  except  in  the  region  near  the  lead¬ 
ing  edge  and  the  apex  where  the  vortex  is  very  close  to  the 
wing  surface  and  where  the  induced  velocities  may  be  expected 
to  be  inaccurate  because  of  the  numerical  problems. 
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In  the  case  of  the  low  aspect  ratio  swept  wing  cer¬ 
tain  guesses  about  the  path  of  the  leading  edge  vortex  can  be 
made  based  on  the  doublet  distribution  on  the  wing  as  described 
in  Section  4.7  of  Reference  4.  However,  a  better  estimate  of 
the  geometry  can  be  obtained  based  on  the  flow  visualization 
studies  and  experimental  pressure  data. 

In  the  case  of  delta  type  wings  with  very  low  aspect 
ratios,  the  geometry  can  be  derived  from  Smith's  conical 
theory  or  experimental  data  if  available.  However,  further 
investigations  into  the  vortex  motion  and  development  of  dif¬ 
ferent  criteria  to  obtain  a  stable  final  geometry  starting 
from  an  arbitrary  initial  location  is  necessary. 

Finally,  it  is  noted  that  execution  of  the  program 
is  extremely  efficient.  For  example,  the  run  time  for  a  low 
aspect  ratio  swept  wing  case  with  three  free  vortices,  192 
boxes  on  the  wing  is  approximately  24  seconds  on  CDC7600 
Cyber  74  system.  A  similar  case  on  the  NASA  Langley  Research 
Center  cluster  requires  44  seconds  on  Cyber  175  series  machine. 

Based  on  the  results  presented  above,  the  following 
conclusions  can  be  made: 

1.  The  program  provides  a  good  correlation  between 
the  predicted  and  observed  total  lift  and  pressure  distribu¬ 
tions  on  the  F-4  type  wings  when  certain  empirical  parameters 
are  known. 

2.  The  progreun  also  provides  reasonably  good  results 
for  delta  type  wings,  but,  certain  refinementr  in  the  prograim 
such  as  inclusion  of  a  leading  edge  sheet  may  be  necessary  to 
improve  the  correlation. 

3.  The  separated  flow  arising  from  the  turbulent  tran¬ 
sition  poses  a  formidable  problem  in  the  cuialysis  within 
the  framework  of  the  potential  theory.  It  is  believed  that 
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certain  problems  introduced  by  the  separated  flow  as  treated 
in  the  analysis  can  be  removed  if  the  determination  of  the 
separated  region  and  the  pressure  there  is  made  after  the 
full  potential  flow  solution  is  obtained.  Such  a  treatment 
of  the  separated  flow  will  include  the  use  of  the  boundary 
layer  theory  and  impressed  pressure  distribution  calculated 
during  the  fully  potential  flow  analysis. 

4.  The  problem  of  the  theoretical  prediction  of 
the  geometry  of  the  leading  edge  vortex  needs  further 
investigation. 

5.  The  program  execution  is  very  efficient  and 
therefore  it  can  be  used  as  a  good  tool  to  obtain  reasonably 
good  results  rather  inexpensively. 


3.  CONCLUSIONS 


Experimental 

On  the  basis  of  the  experimental  investigations  that  were 
conducted  using  wind  tunnel  models,  the  following  general  con¬ 
clusions  Ccin  be  made: 

•  Through  the  use  of  the  controlled  vortex  flows 
interacting  with  a  low  aspect  ratio  swept  lifting 
surface,  a  very  large  increase  in  Cj^  can  be 
achieved. 

•  The  increase  in  C^^  is  accompanied  by  a  destaibi- 
lizing  moment  which,  it  was  shown,  can  be  con¬ 
trolled  with  a  nominal  size  tail  surface  properly 
located  in  the  wake . 

In  addition  to  the  above  general  conclusions,  more 
specific  conclusions  can  also  be  made.  They  are  listed  below. 

1.  With  an  appropriate  combination  of  leading  edge  devices 
such  as  strakes,  snags,  etc.,  the  maximum  lift  coefficient  of 

a  lifting  surface  could  be  increased  by  eibout  60%  as  compared 
to  the  basic  wing  without  these  devices. 

2.  The  lift  curve  slope  for  the  configuration  with  the 

leadingedge  devices  was  linear  with  angle  of  attack  up  to  an 

angle  of  about  29®  whereas  the  lift  curve  slope  for  the  basic 

wing  started  reducing  at  about  15®  reaching  a  minimvim  at  about 

25®  where  occurs, 

max 

3.  Lift  to  drag  ratios  of  the  wing  with  most  advantageous 
vortex  control  devices  did  not  change  substantially  from  that 
of  the  basic  wing  with  no  vortex  control  devices. 

4.  The  most  beneficial  vortex-control  configurations  were 
those  which  involved  at  least  two  interacting  vortices  which 
provided  mutual  stabilization  and  strengthening  as  the  angle 
of  attack  was  increased  above  the  stall  angle  of  the  basic 
wing. 
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5.  The  surface  pressure  distributions  measured  for  the 
various  configurations  indicated  that  the  additional  lift 
obtained  by  the  vortices  interacting  with  the  lifting  surface 
was  primarily  due  to  the  low  pressure  region  in  the  vortex 
cores  acting  on  the  lifting  surface. 

6,.  The  paths  of  the  vortices  over  the  surface  of  the  wing 
were  controlled  primarily  by  the  mutually  induced  effects  of 
the  vortices  and  the  induced  effects  of  the  lifting  surface. 

7.  The  lift  curve  slope  in  the  lower  angle  of  attack 
range  (5®  to  15®)  could  be  increased  significantly  by  inducing 
leading  edge  separation  by  means  of  a  sharp  edge  separator 
plate  along  the  leading  edge. 

8.  External  appendages  having  two  free  edges  which  gener¬ 
ated  vortices  of  opposite  sense  outside  of  the  pressure  field 
of  the  lifting  surface  were  not  effective  in  enhancing  the  wing 
vortex  field  and  thus  the  lift  at  high  angles  of  attack. 

9.  Planform  tailoring  in  the  leading  edge  region  such  that 
only  a  single  free  vortex  is  formed  within  the  wing  pressure 
field  had  the  greatest  potential  for  increasing  the  lift  capa¬ 
bilities  of  a  wing  at  high  angles  of  attack  due  to  favorable 
interacting  vortex  flows. 

10.  Over  a  wide  range  of  wing  angle  of  attack  0®<a^<30®, 

the  tail  effectiveness/  g^/  remained  relatively  constant. 

t 

11.  The  aerodynamic  angle  of  attack  at  the  tail  surface 
never  exceeded  10®  even  though  geometric  angle  with  respect 
to  the  free  stream  was  as  high  as  24®.  This  is  mainly  due  to 
the  induced  effects  of  the  wing  surface  and  vortices  at  the 
tail  surface. 

12.  Depending  on  the  longitudinal  and  vertical  location  of 
the  tail  surface  relative  to  the  wing,  the  tail  drag  can  have  a 
significant  effect  on  the  trim  of  the  total  configuration. 
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Theoretical 


On  the  basis  of  the  results  obtained  using  the  theoretical 
prediction  technique  that  was  developed,  the  following  conclu¬ 
sions  were  drawn: 

1.  In  terms  of  the  theoretical  analysis,  it  was  concluded 
early  in  the  program  that  existing  theories  and  predictive  tech¬ 
niques  were  inadequate  in  predicting  the  performance  character¬ 
istics  of  the  lifting  surfaces  with  leading  edge  vortex  control 
devices . 

2.  Interaction  of  a  vortex  on  the  suction  side  produced 
beneficial  increments  in  the  suction  lift;  the  vortices  formed 
within  the  pressure  field  of  the  vortex  had  a  more  beneficial 
effect  as  compared  to  the  free  vortices  generated  by  devices 
upstream. 

3.  Theoretical  analyses  to  be  used  in  the  prediction 

of  the  loads  and  pressures  on  low  aspect  ratio  swept  or  delta 
wings  should  include  the  treatment  of  the  nonlinear  suction 
effects  and  turbulent  transition  separation  of  the  flow  within 
the  framework  of  a  lifting  surface  theory. 

4.  The  application  of  the  prediction  technique  to  a  low 
aspect  ratio  wing  (F-4  Phantom  wing  planform)  has  shown  that 
fairly  good  correlation  can  be  obtained  between  the  predicted 
and  experimental  pressure  distribution  over  the  wing  planform. 

5.  The  application  of  the  theoretical  prediction  technique 
to  low  aspect  ratio  delta  wings  has  shown  that  only  a  fair  cor¬ 
relation  was  obtained  between  predicted  and  experimental  results. 
The  lack  of  good  correlation  is  believed  due  to  the  absence  of 

a  connecting  vortex  sheet  between  the  core  of  the  vortex  and 
the  leading  edge  in  the  lifting  surface. 

6.  The  application  of  the  prediction  technique  has  also 
shown  the  existence  of  slow  divergence  of  the  predicted  leading 
edge  vortex  geometry.  This  may  be  due  to  the  inaccurate 
modeling  of  the  induced  effects  of  the  mixed  flow  (potential 
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cuid  separated  flows)  over  the  wing  planform,  the  empirical 
treatment  of  the  separated  flow  over  the  wing  or  the  absence 
of  a  roll-up  vortex  sheet  between  the  vortex  core  and  the 
leading  edge. 

7.  In  spite  of  some  limitations  in  the  prediction  tech¬ 
nique,  it  can  be  used  to  generate  realistic  pressure  distribu¬ 
tions  over  low  aspect  ratio  wings  quickly  and  inexpensively 
and  thus  is  useful  for  design  studies  of  advanced  wing 
configurations . 
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4 .  RECOMMEND AT IONS 


Based  on  the  research  investigations  that  were  con¬ 
ducted,  the  following  reconunendations  are  made: 

1.  The  theoretical  analysis  be  modified  to  account  for 
yawed  flight. 

2.  The  theoretical  analysis  be  expanded  to  include  con¬ 
trol  surfaces,  such  as  the  horizontal  and/or  vertical  tail, 
and  the  fuselage  in  order  to  study  the  effect  of  the  strong 
vortex  flows  on  the  static  stability  characteristics  of  air¬ 
craft  configurations. 

3.  The  basic  potential  flow  theory  be  modified  to  in¬ 
clude  the  effects  of  the  leading  edge  vortex  sheet  which 
rolls  up  into  the  leading  edge  vortex. 

4.  An  improved  model  of  the  vortex  core  breakdown  be 
developed  and  incorporated  in  order  to  make  the  analysis  less 
empirical. 

5.  The  vortex  position  iteration  be  altered  in  order 
to  provide  convergence  of  the  vortex  geometry  starting  from 
an  arbitrary  position. 

6.  The  treatment  of  the  separated  flow  be  improved  by 
applying  more  realistic  boundary  layer  separation  criteria. 
This  should  remove  some  of  the  numerical  problems  associated 
with  the  calculation  of  the  separated  flow  regions. 
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Static  pressure 


Chordwise  pressure  distribution  on  the  upper  surface 
of  the  F-4  wing  at  M^*0.7  due  to  a  snag  vortex 
wing  interaction  (from  ref.  6.). 
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FIG.  3.  Simplified  wing  vortex  interaction  model. 
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FIG.  4.  Vortex  lift  as  a  function  of  angle  of  attack  for  the 
strake  wing  configuration. 
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FIG.  8.  Comparison  of  the  incremental  drag  performance 
of  different  high  lift  devices. 
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FIG.  10.  Model  Planform 
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FIG.  18.  Surface  Flow  for  Wing-Strake  Configuration 
a  =  13.1  Degrees. 
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FIG.  22.  Surface  Plow  for  Wing-Strake  Configuration 
a  =  27.7  Degrees 


FIG.  23.  A  Typical  Multi-strake  Configuration. 


FIG.  24.  A  Configuration  with  a  Single  Strake  and  a 
Separator  Plate. 
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PIG.  43.  The  Vortex  Core  Geometry  and  Models. 


FIG.  45.  The  Wing  and  Vortex  Geometry  Configuration  for 
a  ■  21.6  Degrees. 


112 


\ 

\ 

\ 


FIG.  47.  The  Wing  and  Vortex  Geometry  Configuration  for 
a  «  27.7  Degrees. 
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FIG.  49.  Comparison  of  Calculated  Total  Lift  Coefficients  with 
Observed  Values. 
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M’PKNDIX  A 


DKSCKU'TION  OF  THK  COMl'UTKK  PROORAM 


Tho  flow  chart  lor  Lho  pro^jram  soqiwonts  is  prv>vivltHl  in  Fio- 
urt*  A- 1 .  Thoso  soqmonts  are  briofly  doscr  ibod  to  rolato  thorn 
to  tho  appropriato  subroutines  in  tho  program.  A  syml.H'»l  di- 
roctory  is  then  provided  for  tho  computer  variables  rolat inu 
them,  wherever  possible,  to  the  correspotuiii^u  aUiobraic 
variables . 

Thu  input  to  the  pro>.u'Am  is  read  from  cards  in  Program  VOWUA. 

The  calculation  of  the  winij  box  arrangement  and  cot^trol  pi'ints 
is  done  in  svibrovit inos  hOXES  and  CONTROL. 

The  calculations  involved  in  tho  attached  flow  doublet  distri¬ 
bution  solution  are  performed  in  subroutines  AICWINO,  V'F.HNO 
and  SOLVE. 

Tho  pressures  in  the  attached  and  seviaratod  revjiv'ns  are  comv'ut  ed 
in  subroutine  POTENL. 

The  displacement  of  the  free  vortices  is  calculatevl  in  subrout  iiu' 
FREEP. 

Tho  incremental  suction  pressures  due  to  the  free  vortices  are 
calculated  in  subroutine  V’TXPKS.  The  pressure  component  s  at  e 
added  accord  inq  to  F.q .  2. 

Each  of  these  subroutines  calls  certain  other  secondary  routines. 
All  the  subroutines  are  described  in  a  chrono  lovt  i  ca  I  \'>rder  in 
which  they  are  accessed  in  the  proqram. 


SYMBOL  DIRECTORY 


Computer  Symbol 

Algebraic  Symbol 

Description 

AA 

“o 

Geometric  angle  of  attack 

ABEVL 

Bevel  angle  for  delta  wing 
(used  only  if  NDELTA/0) 

ADELTA 

tt/2  -  semiapex  angle  of  the 
delta  wing  (used  only  if 
NDELTA/0) 

AS 

“s 

Empirical  angle  used  in 
the  determination  of  the 
separated  region 

B 

b 

Semispan  of  the  wing 

BSML 

Nondimensional  tangential 
velocity  in  the  spanwise 
direction 

CBAR 

Mean  aerodynamic  chord. 

CHD 

C 

Chord,  local 

CLSAVE 

'  ^L 

W  V 

Lift  coefficient  due  to 
wing  and  suction  resp 

CLTOT 

Total  lift  coefficient 

CMXT 

Moment  coefficients  around 

X  axis  due  to  forces  on  the 
semi -wing 

C^^YT 

^m„'^m  '^m 

y  y„  yv 

Moment  coefficients  around 
y  axis  due  to  forces  on  the 
semi -wing 

CPK 

c 

^att^ 

Pressure  coefficient  at  a 
control  point  on  a  box  in 
the  attached  flow  region 

CPS 

C 

Psep^ 

Pressure  coefficient  at  a 
control  point  on  a  box  in 
the  separated  flow  region 

CPSL 

(dz  /dx) 

i 

Chordwise  slope  of  the  sur¬ 
face  at  a  control  point 

CPTS 

"S 

Suction  pressure  coefficient 
due  to  free  vortices  at  a 
point  on  the  wing 

Computer  Syinbol 
CPXW 


Algebraic  Symbol 


Description 

X  coordinate  of  a  control 
point  on  a  box  on  the  wing 
(ft) 


CPYW 

y  coordinate  of  a  control 
point  on  a  box  on  the 
wing  (ft) 

CPZW 

z 

‘^i 

z  coordinate  of  a  control 
point  on  a  box  on  the 
wing  (ft) 

CORER 

core  radius  at  a  control 
point  on  the  free  vortex 

DELL 

As 

The  length  of  an  element 
of  a  free  vortex  (ft) 

DL1,DL3 

DM1, DM3 

DN1,DN3 

^1'  ^3 

*^1 '  ***3 

"l'  *^3 

Direction  cosines  used  in 
calculation  of  the  induced 
velocity  due  to  a  vortex 
element 

DTIME 

At 

Time  interval  used  in  the 
displacement  of  the  control 
points  on  a  free  vortex  (sec) 

ERR 

c 

Sum  of  squares  of  the  errors 
in  the  coordinates  in  the 
vortex  displacement  calcula¬ 
tions 

ERRP 

Input  parameter  used  in  the 
comparison  for  convergence 
of  vortex  geometry 

FPI 

4it 

FST 

U^Sin  a 
o 

s 

Maximum  allowable  outward 
normal  velocity  at  a  control 
point  on  the  wing  before  the 
flow  on  the  box  containing 
the  control  point  is  considered 
separated  (ft/sec) 

r  Circulation  of  the  free  vor¬ 

tex  elements  in  the  constant 
strength  region  (ft®/sec) 


GAMC 


Description 


computer  Symbol  Algebraic  Symbol 

GAMK  d . 

D 

GAMV  r 

V 

GRC  7D,Vd^ 

GRS 

ISUM 

IT 

I  TRY 

IW, I WRITE 

KINF  r/2  TT 

NACA 


NBL 

^1 

NBU 

NCL 

^c 

^1 

NCP 

^c 

p 

NCU 

Doublet  density  on  a  box  in 
the  attached  flow  region  of 
the  wing  (calculated  by 
satisfying  the  normal  flow 
boundary  condition)  (ft^/sec) 

Circulation  of  the  free  vor¬ 
tex  elements  in  the  variable 
strength  (formation)  region 
(f t^/sec) 

Gradient  of  the  doublet  dis¬ 
tribution  in  the  chordwise 
tangential  direction  (ft/sec) 

Gradient  of  the  doublet  dis¬ 
tribution  in  the  spanwise 
tangential  direction  (ft/sec) 

Parameter  used  to  compute  the 
pressure  difference  between 
the  upper  and  lower  surface 

Iteration  parameter  in  the 
iteration  for  the  doublet 
strength  solution 

Iteration  parcimeter  in  the 
vortex  geometry  iterations 

Output  control  parameters 


Circulation  of  a  vortex 
element  divided  by  2it, 

(ft^/sec) 

Control  parameter  that  effects 
the  computation  of  z,  dz/dx 
according  to  NACA  4  digit  air¬ 
foil  geometry  equation 

Number  of  boxes  on  the  lower 
surface 

Number  of  boxes  on  the  upper 
surface 

Number  of  chordwise  strips  of 
boxes  on  the  lower  surface 

Number  of  chordwise  lines  over 
which  the  incremental  suction 
pressure  due  to  free  vortices 
is  computed 

Number  of  chordwise  strips  of 
boxes  on  the  upper  surface  of 
the  wing 


Computer  Symbol  Algebraic  Symbol 


Description 


NDELTA 

Integer  parameter  to  select 
the  layout  of  points  at  which 
the  incremental  suction  pres¬ 
sures  are  computed 

NEGC 

Number  of  elements  on  each 
free  vortex  in  the  formation 
region 

NINIT 

Number  of  iterations  on  vortex 
geometry  for  each  iteration  on 
doublet  strength 

NIT 

Number  of  iterations  on  the 
doublet  strength  solution 

NKVRT 

The  identification  number  of 
the  vortex  whose  variable 
strength  is  derived  from  the 
doublet  solution  (subroutine 
VSTRNG) 

NOVE 

"el 

Number  of  elements  on  each 
free  vortex  (including  the 
ones  in  the  formation  region) 

NRL 

Number  of  boxes  in  each  chord 
wise  strip  of  boxes  on  the 
lower  surface  of  the  wing 

NRP 

N 

Number  of  spanwise  lines 
(defined  by  x  «  constant  or 
x  =  constant)  over  which  AC 
is  computed  ^v^ 

NRU 

Number  of  boxes  in  each  chord- 
wise  strip  of  boxes  on  the 
upper  surface  of  the  wing 

NSEP 

"‘e 

Number  of  boxes  on  the  trailing 
edge  of  the  wing 

NSOLVE 

"a 

Number  of  boxes  in  the  attached 
region  of  the  wing 

NTOT 

N 

Total  number  of  boxes  on  the 
wing  surface 

NVARV 

Number  of  variable  strength 
elements  <10 

NVORT 

Number  of  free  vortices 

NWBOX 

NTE(i) 

Identification  numbers  of 
boxes  on  the  trailing  edge 

QNET 

The  magnitude  of  the  total 

velocity  at  a  control  point 
on  the  free  vortex  (ft/sec) 
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Computer  Symbol 

Algebraic  Symbol 

Description 

Self-induced  velocity  (in 
the  direction  perpendicular 
to  the  plane  containing  two 
adjacent  elements)  due  to  an 
adjacent  element  of  a  control 
point  on  the  vortex  (ft/sec) 

R 

Radius  of  curvature  of  the 
vortex  at  a  control  point  (ft) 

RBUR 

^bur 

Radius  of  the  core  just  before 
the  vortex  core  bursts  (ft) 

SIG 

Aerodynamic  influence  coeffi¬ 
cients  at  control  points  i 
due  to  unit  density  doublet 
panel  on  boxes  j  and  j' (ft/sec) 

SLLE 

(dx/dy) 

Slope  of  the  leading  edge 
(assumed  straight) 

SLTE 

(dy/dy)^g 

Slope  of  the  straight  trail¬ 
ing  edge 

SMLA 

“i 

Turbulent  eddy  viscosity  param¬ 
eter  at  a  point  on  the  free 
vortex 

SMLD 

d 

The  distance  of  a  point  on  the 
wing  from  a  vortex  in  the  radial 
direction 

SMLH 

d' 

(d/r^) 

SMLVT 

< 

ft 

a 

o 

Chordwise  tangential  velocity 
at  a  controj.  point  on  the  wing 
(nondimensional) 

TAU,TAUP 

T 

Thickness  ratio 

TEDGE 

The  thickness  of  the  edge  of 
the  delta  wing  in  inches 

THC 

®i 

The  angle  between  the  local 
tangent  in  the  chordwise  di- 

rection  and  the  horizontal 
plane  (radians) 

The  thickness  of  the  flat 
section  of  the  A  wing  in 
inches 
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THICK 


Computer  Symbol  Algebraic  Symbol 


Description 


UZERO 

^o 

Free  stream  velocity  (ft/sec) 

VAX 

''a. 

X 

Axial  velocity  at  the  core 
center  at  a  point  on  the 
vortex  (ft/sec) 

VEOX 

VEDY 

(  ^Inew 

[^v  ) 

► 

(x,  y,  z)  coordinates  of  a 
control  point  after  displace¬ 
ment  of  old  point  (i-1) ,  ft 

\  /new 

VEDZ 

(s) 

'  /  new 

i 

VIS AIR 

/ 

V 

Kinematic  viscosity  of  air 
(ft^/sec) 

VPHIND 

-H 

■e- 

> 

Swirl  velocity  at  a  point  on 
the  wing  to  a  free  vortex 
(ft/sec) 

VEX 

('•‘I-.. 

VEY 

> 

(x,  y,  z)  coordinates  of 
control  points  on  the  free 
vortex,  ft 

VEZ 

N..., 

VSTRFAC 

A  factor  used  in  the  computation 
of  the  strength  of  the  vortex  in 
subroutine  VSTRNG 

WI 

Component  of  normal  velocity 
at  a  control  point  on  the 
wing  due  to  the  free  vortices 
(ft/sec) 

WNET 

•H 

c 

External  normal  velocity  at 
control  points  on  the  wing 
(ft/sec) 

XB 

^b 

X  coordinate  of  the  corners  of 
the  boxes  (ft) 
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Computer  Symbol 

Algebraic  Symbol 

Description 

XCAP 

®i 

The  distance  along  the  vor¬ 
tex  centerline  at  control 
points  (ft) 

X£R 

^r 

X 

The  error  in  the  x  coordinate 
of  a  displaced  control  point  on 
the  free  vortex 

XL£ 

x-, 

le 

X  coordinate  of  the  leading 
edge  of  a  given  chord 

XLEAD 

The  parameter  used  in  subroutine 
VTXPRS  to  decide  which  x=constan; 
lines  are  used  to  compute  suction 
pressure  at  (used  only  if  NDELTA^O) 

XMAC 

The  x  coordinate  of  the  quarter 
chord  point  of  the  mean  aero¬ 
dynamic  chord 

XPER 

X 

Nondiroensional  coordinates, 

(given  by  (x-x^^^)/ (x^^-x^^)  , 

giving  the  corners  of  the 
boxes 

XTE 

^te 

X  coordinate  of  the  trailing 
edge  of  a  given  chord  (ft) 

XWP 

X 

X  coordinates  of  the  corners 
of  the  semi -wing  (ft) 

YB 

^b 

y  coordinates  of  corners  of 
boxes  (ft) 

YER 

Error  in  the  y  coordinate 
of  a  displaced  control  point 
on  a  free  vortex 

YMAC 

The  y  coordinate  of  the  quarter 
chord  point  of  the  mean  aero¬ 
dynamic  chord 

YPER 

y 

Nondimensional  y  coordinates 
of  the  corners  of  boxes; 
y  =  y/semispan 

YWP 

^Wi 

y  coordinates  of  corner  of 
the  semi -wing  (ft) 

ZB 

"b 

2  coordinates  of  corners  of 
the  boxes  (ft) 

Computer  Symbol 


Algebraic  Symbol 


Description 


ZER 


ZZERO 


in  the  z  coordinate 
of  a  displaced  control  point 
on  a  free  vortex 

The  transition  distcuice  for 
the  vortex  where  Batchelor's 
model  becomes  applicable  (ft) 


DESCRIPTION  OF  ROUTINES 

PROGRAM  VOWHA  (INPUT,  OUTPUT,  TAPEl,  TAPES  »  INPUT,  TAPES  * 
OUTPUT) 

This  is  the  controlling  segment  of  the  progreun.  It  per¬ 
forms  the  following  functions  by  calling  the  appropriate  sub¬ 
routines  : 

I.  Defines  the  following  files  on  the  scope  operating 
system. 

A.  INPUT  -  the  standard  system  input. 

B.  OUTPUT  -  the  standard  system  output. 

C.  TAPEl  -  unformatted  file  which  stores  the  influence 

coefficients  of  the  doublet  panels  on 
the  wing 

II.  Reads  input  in  the  following  manner: 

UZERO,  AA,  AS,  TAU,  VISAIR,  VLIMIT 
XWP,  YWP 

CBAR,  SMAC,  YMAC 

r^lARB,  IWRITE,  IW,  NDELTA,  NGRAD,  NSIG 

(If  NDELTA  f-  0)  THICK,  TEDGE ,  ADELTA,  ABEVEL 


(If  MARB  f  0)  NRU, 
NSEP,  XB,  YB,  NWBOX 

NCU, 

NRL, 

NCL,  NTOT,  NBU,  NBL, 

(If  ^ARB  «  0)  NRU, 

NCU, 

XPER, 

YPER 

NVORT,  NIT,  NINIT, 

NRP, 

NCP, 

ISUM 

(If  NVORT  f  0)  NOVE,  NEGC,  VEX,  VEY,  VE2 , 
CORER,  (If  NEGC  f  0)  GAMV,  GAMC,  2ZERO 


ERRP,  VILMT,  XLEAD 
VSTRFAC 

NSTR,  NKVRT,  NVARV 
NACA 

III.  Writes  out  input  parameters. 
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IV.  Trtfc  following  subroutines  are  called  in  the  sequential 

order  given  below. 

A.  If  MARB  =  0,  calls  BOXES  to  compute  the  wing  box 
arrangement . 

B.  Calls  CONTROL  to  compute  the  control  points  and 
the  z  coordinates  of  box  corners. 

C.  Calls  AICWING  to  compute  the  aerodynamic  influence 
coefficients . 

D.  Calls  VELIND  to  compute  the  normal  velocity  in¬ 
duced  by  the  free  vortices. 

E.  Calls  SOLVE  to  determine  the  attached  flow  region, 
set  up  the  boundary  conditions  in  the  attached 
region,  and  solve  for  wing  doublet  strengths  d^. 

F.  Calls  POTENL  to  compute  the  pressure  coefficients 
due  to  doublet  panels  in  the  attached  and  separated 
regions  of  the  wing. 

G.  Calls  VSTRNG  to  derive  the  vortex  strength  from  the 
doublet  solution. 

H.  If  NVORT  0  calls  FREEP  to  compute  the  displaced 
vortex  geometry  which  satisfies  the  force  free 
condition. 

I.  If  NIT  >  0  and  IT  <  NIT)  reruns  to  step  C. 

J.  If  IT  >  NIT  and  NVORT  0,  calls  VTXPRS  to  com¬ 

pute  the  incremental  suction  pressures  and  result¬ 
ing  lift. 

X.  Adds  the  force  and  moment  coefficients  due  to 
suction  amd  wing. 
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L 


Subroutine  BOXES 

I.  Computes  the  (x,y)  coordinates  of  the  corners  of 
each  box  for  a  given  set  of  constant  percentage 
chord  (x)  lines  and  constant  percentage  semispan (y) 
lines  into  which  the  semiwing  planform  is  divided. 

These  are  given  by: 

“b  =  -le  ^  -  =‘te> 

Yb  =  yt 

where,  x^^^  and  x^^  are  obtained  for  a  given  y  via 
calls  to  subroutine  EDGES. 

II.  Computes  the  total  number  of  boxes  on  the  upper  sur¬ 
face  (N^) ,  the  total  number  of  boxes  on  the  lower 
surface  (N^)  and  the  total  number  of  boxes  on  the 
semiwing  planform  (N) .  These  are  given  by: 


III.  Computes  the  total  number  of  boxes  at  the  trailing 

edge  and  identifies  the  number  associated  with 

each  trailing  edge  box  (NTE(i)). 

is  given  by; 

te  c^_j 

Subroutine  EDGES 

I.  Computes  the  leading  edge  and  trailing  edge  x 

Xj^g  and  x^^  respectively,  for  a  given  value  of  y. 

are  given  by; 

and 

*te  =  "^^4  (^/dy)te  ' 

Subroutine  CONTROL 

I.  Calls  subroutine  CNTRD  to  compute  the  coordinates 
(Xc»  y^)  of  the  control  point  of  each  box. 

II.  Calls  subroutine  EDGES  to  compute  the  leading  and 

trailing  edge  x  coordinates  and  respectively) 

for  a  given  value  of  y 

c  • 

III.  Calls  subroutine  ZCOR  emd  ZSLP  to  compute  the 

coordinate  and  the  chordwise  slope  (dz/dx)  at  each 
control  point  (x_,  y_) . 

W  w 

IV.  Calls  subroutines  EDGES  and  ZCOR  to  determine  the 
Zjjj  coordinates  of  the  corners  of  all  the  boxes. 


coordinates , 
These 
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Subroutine  CNTRD 


I 


4 


L 


Computes  the  x  and  y  coordinates  of  the  centroid  of 
a  given  box  on  the  wing.  They  are  given  by: 


(^a  • 

^124 

’‘b 

^2  34^ 

/ 

^^124  ^234 

<ya  • 

^124 

Vb  ‘ 

^234^ 

/ 

(Ai24  +  A234 

^124 

and  A234 

are 

areas 

of 

triangles  ( 

(234)  respectively,  and 

Xg  »  (Xj^+X2+X4)/3.  y^  =  (y3^+y2+y3)/3. 
xjj  -  (x2+x4+x3)/3.  =  (y2+y3+y4)/3. 


=  (Zj^+Z2+Z3)/3. 


^  (Z2+Z3+Z4)/3, 


1  (Xj^,yj^,Zj^) 


2  (x2 »y2 » Z2 ) 
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Subroutine  ZCOR 


I.  For  a  given  airfoil  shape  given  by  z  =  t/2  f  (x) , 
computes  the  z  coordinate  at  any  given  point  on 
the  surface  of  the  wing.  It  is  assumed  that  the 
airfoil  shape  doesn't  change  along  the  span.  For  a 
flat  plate  airfoil  of  thickness  t  with  semicircular 
leading  euid  trailing  edges.  The  airfoil  shape  is 
described  by 

f  (x)  «  0,  X  <  0 


f  (X) 
f  (X) 

f  (x) 

f  (X) 


V  (t/2)  ^-(x-t/2)S 

t/2  , 

V  (t/2) ^-(x-1+t/2) ^ 

0, 


0  <  X  <  t/2 

t/2  <  X  <  (I-t/2) 

(I-t/2)  <  X  <  1 
X  >  1. 


Subroutine  ZSLP 

I.  Computes  the  chordwise  slope  (dz/dx)  of  the  surface 
at  any  point  on  the  wing  from 


dz 


t/2  df (x) 
dx 


The  same  shape  of  f(x)  given  in  subroutine  ZCOR  is  used 

to  compute  .  The  in  the  various  ranges  are 

dx  dx 


given  below: 

df (x)  ^ 

dx  V  (t/2)  * 


-  X  _ 

(X  -  t/2) ^ 


0  <  X  <  t/2 


0 


t/2  <  X  <  1  -  t/2 


,  -  (X  -  1  /2) _ 

V(t/2)2  -  (X  -  1  +  t/2^ 


(I-t/2)  <  X  <  1 


»  0 


X  >  1. 
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Subroutine  MOWING 


I.  Computes  the  aerodynamic  influence  coefficients, 
o^j,  due  to  the  box  j  and  its  image  j'  at  the 

control  point  of  box  i  (see  the  figure  below) . 


Calls  subroutine  UVW  to  compute,  ^ij' 

cartesian  velocity  components  at  point  i  due  to  boxes 
j  emd  j'.  Confutes  using  the  following  formula 

“ij  ■  “ij  ®i  -  "ij  =°*  ®i 


where  0^  is  given  by 

e. 


tem”^  (dz/dx) 


When  the  box  j  and  j'  are  at  the  trailing  edge  an 
additional  contribution  due  to  trailing  wake  is  also 
included  in  u^^  and  These  are  similarly  cal¬ 

culated  using  subroutine  UVW. 

II.  Stores  the  on  TAPE  lif  the  number  of  iterations 

is  greater  than  1. 


Subroutine  UVW 

I.  Computes  the  velocity  components  (u,  v,  w)  due  to 

a  straight  vortex  element  AB  (see  figure)  of  strength 
r  at  a  point  P. 
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where 

AP  and  BP  are  lengths  of  lines  AP  and  BP,  n^^) 
and  m^/  n^)  are  the  direction  cosines  and  0^  is  in¬ 
cluded  angle  given  by 


Subroutine  VELIND 

I.  Computes  the  Cartesian  components  (u^,  v^,  w^)  at 

control  point  i  due  to  the  free  vortices  whose 
geometry  and  strength  is  known  by  calling  subroutine 
UVW. 


II.  Computes  the  normal  component  of  velocity  due  to  free 
vortices  from 

w„  =  u.  Sin  0.  -  w.  Cos  0. 
n^^  1  1  1  i 
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Subroutine  SOLVE 

I.  Computes  the  net  normal  velocity  component,  ,  at 

a  control  point  i,  due  to  free  stream  velocity  and 
free  vortices.  It  is  given  by: 
w.,  =  Sin  (a^  -6.)  -  w„ 


Determines  the  attached  flow  region  by  conparing  the 

outward  normal  w^  with  Sin  a„.  The  box  i  is 

n^  o  s 

considered  attached  when 


<  Sin 
o  s 


outward 


III.  Sets  up  equations 


E  “ij  =  “n,  ^ 


where  is  the  number  of  boxes  where  flow  is  attached. 

IV.  Solves  the  cUsove  equation  by  calling  subroutine  SLSESO 
for  values  of  dj 

Subroutine  SLSESO 


I.  Solves  the  linear  algebraic  system  of  equation 

i  =  1,  •  •  •  N 

row  by  row. 

The  for  (j  =  1,  N)  are  obtained  from  subroutine 
READR. 


Subroutine  READR 

I.  Reads  the  influence  coefficients  from  TAPE  1. 

II.  If  any  separated  boxes  are  present,  automatically 
eliminates  the  rows  of  <7^^  related  to  those  boxes 
from  computation. 
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Subroutine  VSTRNG 

I.  This  subroutine  interpolates  between  the  values  of 
difference  of  doublet  strengths  on  the  lower  surface 
from  that  on  the  upper  surface  to  obtain  strength  of 
the  leading  edge  vortex.  (This  was  used  for  NDELTA  ^  0) . 

Subroutine  FREEP 

1.  Computes  the  core  radii  at  the  end  points  of  all  free 
vortex  elements  using  the  following  formulas. 

r^  =  .0264  (Sq+s^)  ®  ^  ®i  ^ 

rc^  =  .0264  ^ 

where. 

Sj_  is  the  distance  along  the  leading  edge  vortex,  at 
the  free  vortex  element  i, 

Zq  is  the  value  of  at  which  the  viscous  vortex  core 
model  is  chcmged. 


and  is  the  vortex  strength  of  the  free  vortex  element. 

The  eddy  viscosity  parameters  and  are  obtained 
via  calls  to  subroutine  EDDYV. 

The  axial  velocity  in  the  core  of  the  free  vortex 
element  i  is  given  by: 


A-21 


1.26  r 


0.02487  r 


The  point  of  bursting  of  free  vortex  is  specified  by 
the  location  along  the  vortex  at  which  the  core  radius 
equals  r^^  and  is  given  by 


.04  (  s^  -t-  s^) 


For  vortex  elements  beyond  the  point  of  bursting,  the 
axial  core  velocity  is  given  by 

.002488  r..  a..  /s_  +  z_ 


II.  Calls  subroutine  VELBOX  to  confute  the  induced  velocity 
components  due  to  attached  doublet  distribution. 

III.  Computes  the  induced  velocity  due  to  all  the  elements  of 
the  free  vortices  at  a  control  point  i  on  the  vortex 
except  due  to  the  adjacent  elements. 

[V.  Calls  subroutine  SLFUVW  to  compute  the  effect  of  the 
adjacent  elements  a  control  point  on  the  vortex. 

V.  Limits  the  induced  contribution  due  to  the  doublets  and 
the  vortices  to  a  given  value  VILMT  provided  as  input. 

VI.  Confutes  the  total  velocity  by  adding  all  the  velocity 
components  (including  the  free  stream  component)  to 
obtain 


u.  i  v^  j  +  w,  k  . 


1 


Computes  the  time  interval,  iit,  required  for  displace¬ 
ment  of  the  point  i  which  is  given  by 
AS.. 


Ati 


where 


VII, 


and 

|q|  =  V  u^^  +  V. ‘  +  w^‘ 

Computes  the  displacements  of  the  point  i  and  adds 
them  to  its  old  coordinates  to  get  the  new  location  of 
point  i  (i+l  in  the  next  iteration) 

=  lx..  \  +  uAt 


'i+l( 


new 


V. 

1 


old 


^■^^Uew  I  ^lold 


Z  \  a  Z 

^■‘■^Uew  ^)old 


+  vAt 


+  wAt 


VIII.  Computes  the  errors  in  the  locations  of  the  points  on 
the  vortex  as  given  by 

E 

I  v,^ii  I  v,^,i 

old 


old 


^  new  I 

I  new 

h-l 

new  ^ 

old 
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Li 


and  sums  the  squares  of  those  errors  as  follows: 


e 


+  E 


IX. 


where,  is  the  number  of  elements  in  each  vortex. 

The  quantity is  compared  with  (an  input) 

for  convergence.  The  geometry  is  considered  to  have 
converged  if 


<  e  . 
P 


Subroutine  EDDYV 

I.  Computes  the  turbulent  eddy  viscosity  parameter 
from  the  formula 


=  197VRe”  +  .6 


where 

Re 


r. 

2it  V 


where 


V  is  kinematic  viscosity  and  is  the  vortex 

strength  of  the  element  i. 

Subroutine  VELBOX 

I.  Computes  the  velocity  components  (u,  v,  w)  at  a  given 
point  in  space  due  to  the  doublet  distribution  on  the 
attached  region  of  the  wing.  The  computation  proceeds 
essentially  in  the  same  manner  as  in  subroutine  AICWING, 
except  here  the  doublet  strengths  are  known  and  the 
contribution  due  to  all  the  boxes  is  added  together. 

The  image  boxes  and  the  trailing  wake  elements  are 
considered  in  a  similar  manner  as  in  subroutine  AICWING. 
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Sub  routine  SLFUVW 


I.  Computes  the  velocity  components  at  a  point  i  on  the 
free  vortex  due  to  the  adjacent  element  (see  figure 
below) 


The  velocity  perpendicular  to  the  plane  containing  the 
points  (i-1,  i,  i+D  is  given  by 


This  is  the  contribution  from  element  (i,  i+1) .  The 
is  the  radius  of  curvature,  and  (J>^  is  the  included 


II.  Computes  the  appropriate  Cartesian  components  of 


Subroutine  VTXPRS 

I.  Confutes  the  suction  pressures  at  several  points  on 
the  upper  surface  of  the  wing  which  are  selected 
depending  on  value  of  NDELTA. 

If  NDELTA  »  0  the  points  lie  along  lines  defined  by 
X  =  constant.  If  NDELTA  “  1  the  points  lie  on  lines 
defined  by  x  =  constant.  There  are  N  such  lines  and 

P 

on  each  line  there  are  N  points  at  equal  intervals. 


At  each  selected  point,  computes  the  radial  distance, 
d,  of  the  point  from  the  axis  of  the  vortex. 

Computes  the  incremental  suction  pressures  using  the 
following  formulas . 


.672r^ 


-1.32  ri/r 


(13.14065)  (.61)  $ 


A  •  V/  O 


+  .549  /  1  -1.482  e 


0  <  S.  <  Z 
1  o 
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2 


2 


-4 


z 


o 


,02847r^Qi^ 
7i  ~ 


's  +  z  * 

Id.  ,  \* 

o  o 

Y“  2 

(2.48) 

.78  e’*^^ 

^c. 

L  1  J 

\ 

V. 


V. 


Computes  the  total  lift  coefficient,  C  ,  due  to 

- 

V 

incremental  suction  pressures. 

Computes  the  total  moment  coefficients,  C  and  C  , 

'mm' 

about  X  and  y  axes  respectively,  due  to  forces  on  the 
semi-wing. 


VII.  Sets  up  the  arrays  for  interpolation  of  potential  pres¬ 
sure  by  calling  subroutine  INTARR. 

Interpolates  between  the  potential  pressure  values  on  the 
upper  surface  and  adds  the  potential  and  suction  components 
of  the  pressure. 


Subroutine  INTARR 

I.  Sets  up  the  interpolating  arrays  of  values  associated 

with  functions  such  as  pressure  and  doublet  strengths  at 
control  points  on  the  wing.  These  arrays  are  used  to 
interpolate  for  values  of  these  functions  at  points  on 
the  wing  other  than  the  control  points. 
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Subroutine  PQTENL 

I.  Calls  subroutines  VELBOX  to  compute  the  induced  velocity 
components  (u^,  ) ,  due  to  the  doublet  panels  on 

the  attached  flow  region  of  the  wing. 

II.  Computes  the  chordwise  gradients,  Vd^ ,  of  the  doublet 
strength,  by  calling  subroutine  GRADNT. 

III.  Computes  the  chordwise  and  the  spanwise  tangential 
velocity  components  as  given  by 

■  u^  Cos  +  w^  Sin  6^  +  '7d^/2. 


IV.  Computes  the  pressure  coefficient  at  a  control  point 
on  a  box  in  the  attached  flow  region  by  using  the 
Bernoulli  equation  in  the  form 


V.  Computes  the  pressure  coefficients  on  the  boxes  where 
flew  is  separated  using  the  following  formula 


where  |q^|  is  the  total  tangential  velocity  and  does 
not  include  the  local  contribution  Vd^/2. 
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VI.  Computes  the  total  lift  and  moment  coefficients  , 

^m  '  ^m  forces  on  the  semi -wing. 

Subroutine  OUTPUT 

I.  Divides  the  boxes  in  the  chordwise  and  spanwise  lines 
and  prints  out  the  array  containing  the  values  related 
to  the  control  points  of  boxes  on  the  wing  in  those 
chordwise  and  spanwise  lines,  both  for  upper  and  lower 
surfaces. 

Subroutine  GRADNT 

I.  Computes  the  value  of  the  gradient  of  doublet  strength 
at  a  point  on  the  wing. 
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INPUT  OUTLINE 


First 

Card;  (16A5) 

a. 

80  character 

title 

Next 

Card:  (8F10.6) 

a. 

UZERO, 

the  free  stream  velocity,  ft/sec 

b. 

AA, 

the  angle  of  attack,  radians 

c. 

AS, 

the  angle  used  in  the  determination  of 
separated  flow  region,  radians 

d. 

TAU, 

the  thickness  ratio  of  the  wing 

e. 

VISAIR, 

the  kinematic  viscosity  of  air,  ft^/sec 

Next 

Cards  (8F10.5) 

a. 

XWP, 

the  X  coordinates  of  the  four  corners 
of  the  semi-wing,  ft 

b. 

YWP, 

the  y  coordinates  of  the  four  corners 
of  the  semi -wing,  ft 

Next 

Card:  (8F10.5) 

a. 

CBAR, 

the  mean  aerodynamic  chord 

b. 

XMAC, 

the  X  coordinate  of  the  quarter  chord 
point  of  the  mean  aerodynamic  chord 

c . 

YMAC, 

the  y  coordinate  of  the  quarter  chord 
point  of  the  mean  aerodynamic  chord 

Next 

Card:  (1615) 

a. 

MARB, 

the  integer  parameter  controlling  the 
calculation  of  the  x  and  y  coordinates 
of  the  corners  of  boxes  and  related  input 

■  1  if  box  coordinates  are  to  be  read  in 
directly 

»  0  if  the  box  coordinates  are  to  be  com¬ 
puted  in  the  program 

b. 

IWRITE , 

the  integer  parameter  controlling  the  print 
output  of  subroutines  FREEP  and  INTERP 

c. 

IW, 

the  integer  parameter  controlling  the  print 
output  of  subroutines  BOXES,  CONTROL  and 
VTXPRS 

d. 

N DELTA, 

the  integer  parameter  controlling  the  grid 
points  at  which  the  suction  pressure  is 
calculated 

=  l_if  these  points  are  to  lie  on  constant 
X,  y  lines 

*0  if  these  points  are  to  lie  on  constant 

X  lines 

e . 

NGRAD, 

the  integer  parameter  to  control  the 
computation  of  the  gradient  of  doublet 
strength 

y  0  if  the  gradient  is  to  be  computed 
in  both  spanwise  and  chordwise 
direction 

*  0  if  the  gradient  is  to  be  computed 
in  only  the  chordwise  direction 

f . 

NSIG, 

the  integer  parameter  to  control  the 
computation  of  in  subroutine  AlCWING 

=  1  if  the  file  TAPE  1  already  has  the 
stored 

=  0  if  the  are  to  be  newly  computed 

in  the  present  run 

Next 

Card:  (8F10.5) 

To  be  read  if  NDELTA  f  0 

a. 

THICK, 

thickness  of  the  flat  plate  delta  wing, 
inches 

b. 

TEDGE , 

thickness  at  the  edge  of  the  beveled 
delta  wing,  inches 

c. 

ADELTA, 

the  angle  between  the  y  axis  and  the 
leading  edge  of  the  delta  wing,  degrees 

d. 

ABEVEL, 

the  bevel  angle,  degrees 

Next 

Card:  (1615) 

To  be  read  if  MARB  »  0 

a . 

NRU, 

the  number  of  boxes  in  a  chordwise  strip 
of  boxes  on  the  wing 

b. 

NCU, 

the  number  of  chordwise  strips  of  boxes 
on  one  surface  of  the  wing 

Next 

Card:  (16F5.3) 

To  be  read  if  MARB  =  0 

a. 

XPER, 

the  X  values  that  form  spanwise  edges 
of  the  boxes 
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Next  Card;  (16F5.3)  To  be  read  if  MARB  »  0 


a. 

YPER, 

the  y  values  that  form  the  chordwise 
edges  of  the  boxes 

Next 

Card : 

(1615) 

To  be  read  if  MARB  ft  0 

a. 

NRU, 

the  number  of  boxes  on  each  chordwise 
strip  on  the  upper  surface 

b. 

NCU, 

the  number  of  chordwise  strips  of  boxes 
on  the  upper  surface  of  the  wing 

c. 

NRL, 

the  number  of  boxes  on  each  chordwise 
strip  on  the  lower  surface  of  the  wing 

d. 

NCL, 

the  number  of  chordwise  strips  of  boxes 
on  the  lower  surface  of  the  wing 

e . 

NTOT, 

the  total  number  of  boxes  on  the  wing 
(must  be  less  than  200) 

f. 

NBU, 

the  total  number  of  boxes  on  the  upper 
surface  *  NRU  x  NCU 

g- 

NBL, 

the  total  number  of  boxes  on  the  lower 
surface  =  NRL  x  NCL 

h. 

NSEP, 

the  total  number  of  boxes  at  the  trailing 
edge  of  the  wing 

Next 

NTOT  Cards; 

(8F10.5)  To  be  read  if  MARB  ^  0 

a. 

XB, 

the  X  coordinates  of  the  4  corners  of 
each  box,  ft 

b. 

YB, 

the  y  coordinates  of  the  4  corners  of 
each  box,  ft 

Next 

Card; 

(2014) 

To  be  read  if  MARB  ft  0 

a. 

NWBOX, 

the  identification  numbers  of  the  NSEP 

trailing  edge  boxes 
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Next  Card:  (1615) 


a. 

NVORT , 

the  number  of  vortices 

b . 

NIT, 

the  number  of  iterations  on  the  surface 
doublet  distribution 

c. 

NINIT, 

the  number  of  iterations  on  the  vortex 
geometry  for  each  solution  of  the  doublet 
distribution 

d. 

NRP, 

the  number  of  spanwise  lines  on  which  the 
incremental  suction  pressure  is  calculated 

e . 

NCP, 

the  number  of  points  on  each  spanwise  line 
at  which  the  incremental  suction  pressure 
is  calculated 

f . 

ISUM, 

integer  parameter  that  controls  the  cal¬ 
culation  of  the  pressure  (to  be  used  only 

for  very  thin  wings) 

*  1  if  the  pressure  difference  between 
the  upper  and  lower  si  face  is  desired 

=  0  if  the  pressure  on  each  surface  is 
desired 


The  following  NVORT  sets  of  cards  to  be  read  if  NVORT  /  0 
Next  Card:  (1615) 

a.  NV,  the  number  of  points  on  the  vortex, 

NOVE  »  NV  -  1 

b.  NEGC,  the  number  of  elements  on  the  vortex 

with  variable  strength 

Next  NV  Cards :  8F10.5) 


a. 

VEX, 

the  X  coordinates 
vortex,  ft 

of 

the 

point 

on 

the 

b. 

VEY, 

the  y  coordinates 
vortex,  ft 

of 

the 

point 

on 

the 

c. 

VEZ, 

the  z  coordinates 
vortex,  ft 

of 

the 

point 

on 

the 

d. 

CORER, 

the  core  radius  of 

the 

point  on 

the 

vor 

tex,  ft.  Only  the  first  value  needs  to 
be  given,  the  other  values  are  computed 
within  the  program 
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Next  NEGC  Cards  are  present  only  if  NEGC  p  0:  (8F10.6) 


a. 

GAMV, 

the  strength  of  the  element  of  vortex, 
in  the  variable  strength  region,  ft^/sec 

Next 

Card: 

(8F10.6) 

a . 

GAMC, 

the  strength  of  constant  portion  of  the 
vortex,  ft^/sec 

b. 

ZZERO 

9 

the  transition  distance,  z^,  ft 

Next 

Card: 

(8F10.6) 

To  be  read  after  NVORT  sets  of  cards  are 
read  in  and  if  NVORT  ^  0 

a. 

ERRP, 

the  error  parameter  for  comparison  in 
subroutine  FREER 

b. 

VILMT 

t 

the  limit  on  induced  velocity 

c. 

XLEAD 

9 

is  the  distance  that  controls  the  lines 
at  which  the  computation  of  suction  pres¬ 
sure  is  done,  (takes  effect  only  when 
delta  wing  is  considered) ,  inches 

Next 

Card: 

(8F10.6) 

a. 

VSTRFAC, 

the  factor  which  is  used  to  derive  the 
strength  of  the  vortex  from  the  doublet 
solution  (to  be  used  only  for  delta  wing) 

Next 

Card: 

(1615) 

a.  NSTR,  the  integer  to  control  whether  the  strength 

derived  from  the  doublet  solution  is  used 
in  the  further  calculations 

0  if  yes 

■  0  if  no 

b.  NKVRT,  the  identifying  number  of  the  vortex  for 

which  the  strength  is  derived  from  the 
doublet  strength 

(has  to  be  1,  2,  3  or  NVORT) 

usually  the  leading  edge  vortex 
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c.  NVARV 


the  number  of  elements  in  the  above 
vortex  whose  strength  is  derived  from 
the  doublet  solution 


Mext  Card;  (1615) 

a.  NACA  integer  to  control  the  computation  of  z 

coordinate  and  dz/dx 

*  1  if  NACA  4  digit  aircoil  equations  are 
to  be  used 

=0  if  the  beveled  delta  equations  are 
to  be  used 

The  format  in  which  the  card  is  read  is  placed  in  the  parenthesis 
next  to  the  card  identifier. 


The  following  limitations  are  placed  on  certain  variables. 


NTOT  <  200 
NRUXNCU  <100 
NSEP  <  30 
%  <  10 

due  to  the  dinvensioning  in  the  program. 

The  following  limitations  are  placed  on  the  variables  due  to 
output  formatting 

NCU  <  10. 


The  vortex  geometry  is  supplied  from  the  experimental  investiga¬ 
tions  in  the  literature.  The  strengths  GAMC  and  GAMV  are  pro¬ 
vided  based  on  the  following.  The  program  can  be  run  with  high 
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vortices  to  obtain  a  potential  solution 

without  separation.  The  difference  of  doublet  strengths  of 
the  leading  edge  boxes  on  the  upper  and  lower  surface  can 
be  plotted  against  y. 


See  the  figure.  The  maximum  of  can  be  used  as  the 

initial  guess  for  GAMC  for  the  leading  edge  vortex.  The  values 
of  GAMV  (the  variable  strengths)  can  be  derived  from  the  same 
plot  based  on  y  coordinates  of  the  leading  edge  vortex.  The 
values  of  these  strengths  for  strake  and  tip  vortices  are  arrived 
at  from  experience. 
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OUTPUT  OUTLINE 


I .  Input 


Computer  Name 

Algebraic  Name  Description 

UZERO 

«o 

free  stream  speed 

AA 

“o 

euigle  of  attack 

AS 

“‘s 

angle  \ised  in  the  determination 
of  the  separated  flow  region 

XWP,YVJP 

*w  '^w 
i  i 

the  X  and  y  coordinates  of  the 
semi-wing  corners 

NVORT 

number  of  vortices 

NIT 

number  of  iterations  on  the 
doublet  distribution  solution 

NINIT 

number  of  iteration  on  the  vor¬ 
tex  geometry 

II.  The  wing 

box  arr2mgement 

information  (MARB  »  0) 

NTOT 

N 

number  of  boxes  on  the  wing 

NBU 

Nu 

number  of  boxes  on  the  upper 
surface 

NBL 

**1 

number  of  boxes  on  the  lower 
surface 

NSEP 

^te 

number  of  boxes  at  the  trail¬ 
ing  edge 

NWBOX 

NTE(i) 

numbers  identifying  the  boxes 
at  the  trailing  edge 

NRU 

number  of  boxes  on  each  chord 
wise  strip  of  boxes  on  the 
upper  surface 

NCU 

number  of  chordwise  strips  of 

boxes  on  the  upper  surface 

number  of  boxes  in  each  chord 
wise  strip  of  boxes  on  the 
lower  surface 


NRL 


Computer  Name 
NCL 


Algebraic  Name 


XB ,  YB »  ZB 

CP2W,CPYW, 

CPZW,CPSL 


Description 


number  of  chordwise  strips 
of  boxes  on  the  lower  surface 


the  coordinates  of  corners 
of  the  boxes 


the  coordinates  and  slopes 
at  the  control  points 


III.  Attached  flow  doublet  distribution  solution 


IT 


the  iteration  number 


GAMK 


tij  the  values  of  uniform  density 

doublet  panel  strength  on  the 
boxes  where  flow  is  attached 

the  difference  of  doublet 
■’  strengths  on  upper  and  lower 

surface 


IV.  Pressure  calculation  in 
PX,PY,PZ 

WORT 

VNOR 

WI 

"i 


the  attached  and  separated  region 

velocity  components  induced 
by  the  doublets  on  the 
attached  flow  region  at 
control  points 

the  velocity  components  in¬ 
duced  by  the  free  vortices 
at  control  points 

the  normal  velocity  induced 
at  the  control  point  due  to 
doublets 

the  normal  velocity  due  to 
external  source  at  control 
points 
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.a*!.*  iJf 


Computer  Name 


Algebraic  Name 


Description 


VELL 

tangential  velocity  at  the 
control  points  in  chordwise 
direction  induced  by  the 
doublets 

WLL 

tangential  velocity  at  the 
control  points  in  chordwise 
direction  induced  by  the 
free  vortices 

UZC 

component  of  the  free  stream 
in  the  chordwise  tangential 
direction  at  the  control 
points 

GRC 

VD 

local  gradient  of  the  doioblet 
distribution 

MP 

number  of  the  control  point 

XP , YP , ZP 

x,y,z 

the  nondimensional  coordinates 
of  the  control  point 

VCHORD 

induced  component  of  velocity 
in  the  chordwise  tangential 
direction 

VSPAN 

''t 

s 

induced  component  of  the 
velocity  in  the  spanwise 
tangential  direction 

CPK 

Satt 

the  pressure  coefficient 
at  the  box  where  the  flow 
is  attached 

CPS 

^P 

sep 

the  pressure  coefficient  on 
the  box  where  flow  is 
separated 

ABOX 

representative  area  of  the 
box 

GAMK 

the  doublet  panel  strength 

CLPU 

the  lift  coefficient  due  to 
upper  surface  pressures  on 
attached  region 

CLSU 

the  lift  coefficient  due  to 

upper  surface  pressures  on 
the  separated  region 
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CMXU,CMyU 

CLSAVE 

CLSTAL 

CMX,CMY 

CLTOT 

CMXT,  CMYT 
XCOP,YCOP 

Vortex  displacement  information 
UU,  W,  WW 

VEX,  VEY,  VEZ 

VEDX,VEDY,VEDZ 

XER, YER, ZER 

CORER  r 

°i 

GAMC,  GAMV 


moment  coefficients  around 
X  and  y  axes  due  to  pres¬ 
sures  on  the  upper  surface 
of  the  semi -wing 

lift  coefficient  due  to  pres¬ 
sures  on  both  surfaces  of 
the  wing  in  the  attached 
region 

lift  coefficients  due  to 
pressures  on  both  surfaces 
of  the  wing  in  the  separated 
region 

moment  coefficients  due  to 
pressures  on  both  the  sxir- 
faces  of  the  wing 

total  lift  coefficient 

total  moment  coefficients 
around  x  and  y  axes 

center  of  pressure  coordi¬ 
nates  for  pressures  acting 
on  the  semi -wing 


the  components  of  the  velocity 
at  control  points  on  the  vortex 
due  to  wing  doublets,  free 
stream,  vortices  themselves, 
auid  axial  velocity  in  the  core 

(x,  y,  z)  coordinates  of  the 
control  points  on  the  vortex 
(old  position) 

(x,  y,  2)  coordinates  of  the 
control  points  on  the  vortex 
(new  position) 

errors  in  the  coordinates  of 
the  position  on  the  vortex 

core  radius  of  the  vortex  at 
the  control  point  on  the  vortex 

the  strength  of  the  elements 
of  the  vortex 
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UNET,VNET,  VvliET 

XERL,YERL,2ERL 


the  total  velocity  components 
at  the  control  points  on  the 
vortex  due  to  all  sources 

total  cumulative  error  squared 
for  the  given  vortex 


V.  Suction  peaks  information 


XF,YP,ZP 

coordinates  of  the  points 
at  which  the  incremental 

suction  pressures. are 
calculated 

XP1,YP1,ZP1 

nondimensional  x,  y,  z 
coordinates  of  the  above 
points 

OPTS 

AC 

incremental  suction  pres¬ 

*^v 

sure  coefficient 

CPAT 

wing  pressure  coefficient 
(interpolated  value) 

CPTOT 

'"‘'tot 

total  C  value 

P 

NOTE; 


Appropriate  f (x)  and  df (x)  have  to  be  substituted  in 

dx 

subroutines  ZCOR  and  ZSLP  to  match  the  airfoil  of  the  wing  for 
which  the  calculations  are  being  made. 


Presently  the  f(x)  and  df (x) 

dx 


four  digit  NACA  airfoils. 


correspond  to  symmetric 


f(x)  is  given  by 

f(x)  =  10  (.2969V^  -  .126x  -  .  3516x^  +  .2843x*  -  .lOlSx'*) 

df(x)  «  10  (.14845/VT  -  126  -  .7032x+  .8529x^  -  .406x’) 
dx 
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